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ABSTRACT 

Wave  propagation  in  a  random,  inhomogeneous  ocean  is 
treated  as  transmission  thru  a  linear,  time-variant,  space- 
variant,  random  communication  channel.   A  consistent  notation 
(vis-a-vis  ad  hoc),  fundamental  input-output  relations,  and 
various  time-space  transformations  for  both  deterministic  and 
random  linear,  time-variant,  space-variant,  filters  are  esta- 
blished.  Using  the  method  of  separation  of  variables  and  the 
W.K.B.  approximation,  a  time-invariant,  space-variant,  random 
transfer  function  of  the  ocean  volume  is  derived.   The  ocean 
volume  is  characterized  by  a  random  index  of  refraction  which 
is  a  function  of  depth.   The  index  of  refraction  is  decomposed 
into  a  deterministic  component  and  a  zero  mean  random  compo- 
nent.  In  addition,  two  example  calculations  are  made.   The 
first  example  involves  the  derivation  of  the  equations  for 
the  random,  output  electrical  signals  at  each  element  in  a 
receive  planar  array  of  complex  weighted  point  sources  in 
terms  of  the  frequency  spectrum  of  the  transmitted  electrical 
signal,  the  transmit  and  receive  arrays,  and  the  transfer 
function  of  the  ocean  medium.   The  second  example  involves 
the  derivation  of  the  coherence  function,  i.e.,  the  auto- 
correlation function  of  the  transfer  function. 


I.   INTRODUCTION 

Since  the  wave  equation  for  small  amplitude  acoustic  sig- 
nals is  linear ,  we  can  represent  the  ocean  medium  as  a  linear , 
time-variant,  space-variant,  random  filter  (system  or  communi- 
cation channel)  in  general.   With  this  interpretation  in  mind, 
refer  to  Fig.  1  which  illustrates  a  basic  bistatic  communica- 
tion channel  geometry,  and  Fig.  2,  which  is  a  mathematical 
block  diagram  representation  of  Fig.  1.   With  respect  to  Fig.  1, 
both  the  transmit  and  receive  apertures  (arrays)  are,  in  gen- 
eral, volume  apertures  (arrays)  and  in  motion.   Before  pro- 
ceeding further,  a  word  of  caution  concerning  Fig.  2.   Note, 

for  example,  that  X„  £   XD^,  Y^   ±   X„,H„,  and  Y  ^  Y,.D_  in  general. 
r  M      T    M     M  M  MR 

The  equations  required  for  describing  the  filter's  input-output 
relationships  and  for  coupling  the  transmitted  and  received 
electrical  signals  to  the  medium  via  the  transmit  and  receive 
apertures  are  developed  in  Sections  II  and  III,  respectively. 

Let  us  now  describe  the  notation  used  in  Fig.  2.   The  posi- 
tion vectors  r   and  r  refer  to  spatial  coordinates  (x  ,y  ,  z  ) 
— o     —  v  o  -*o  o 

and  (x,y,z),  respectively,  and  t  refers  to  time  in  sec.   The 
parameters  f  and  n  are  frequencies  in  HZ.   where  f  represents 
input  or  transmitted  frequencies  while  n  represents  output  or 
received  frequencies.   Note,  that  if  n  ^   f,  Doppler  spread  is 
implied. 

The  quantities  a,  v_,  3L,  and  y_  are  vectors  whose  components 
are  spatial  frequencies  with  units  of  cycles/m.  Since  spatial 
frequencies  are  related  to  both  direction  cosines  and  wavelength, 
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Y  (DEPTH) 


Fig.  1.   Basic  bistatic  communication  channel 
geometry. 
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and  hence,  wavenumber  components,  they  represent  directions 
of  wave  propagation.   The  vector  v_  represents  input  or  trans- 
mitted spatial  frequencies  into  the  medium  as  in  XM(f,v_), 
while  jJ  represents  output  or  received  spatial  frequencies 
from  the  medium  as  in  Y  (n  ,3_)  .   Note,  that  if  3_  ^  v_f  angular 
spread  (scatter)  is  implied. 

The  remaining  expressions  found  in  Fig.  2  are  further 
described  in  the  following  list: 

x(t,r)  -  input  electrical  signal  to  transmit  electro- 
acoustic  transducer  applied  at  time  t  and 
spatial  location  r  of  transducer. 
X(f,o_)  -  frequency  (f)  and  angular  (a_)  spectrum  of  input 

electrical  signal. 
A_(f,r)-  complex  frequency  response  at  spatial  location 
r  of  transmit  transducer.   Also  referred  to  as 
the  complex  transmit  aperture. 
DT(f,oO-  transmit  far-field  directivity  function  or  beam 

pattern. 
xM(t,r)-  input  acoustic  signal  to  the  medium  applied  at 
time  t  and  spatial  location  r.   Also,  output 
acoustic  signal  from  transmit  electro-acoustic 
transducer. 
X-.(f,v_)-  frequency  (f)  and  angular  (_v)  spectrum  of  input 

acoustic  signal, 
h  (x,r  ;  t,r)  -  time-variant,  space-variant  impulse  re- 
sponse of  the  ocean  medium.   It  represents  the 


response  of  the  medium  at  time  t  and  spatial 
location  r  due  to  the  application  of  an  unit 
impulse  at  time  Ct-T)sec,,  or  t  sec.  ago,  at 
a  distance  |r-r  |  m.  away  (see  Fig.  1). 

H  (f,_v;t,r)  -  time-variant,  space-variant  transfer 
function  of  the  ocean  medium. 

y  (t,r)  -  output  acoustic  signal  from  the  medium  at  time 
t  and  spatial  location  r.   Also,  input  acoustic 
signal  to  receive  electro-acoustic  transducer. 

Y  (n,£)  -  frequency  (n)  and  angular  (B_)  spectrum  of  out- 
put acoustic  signal. 

A-,(ri,r)  -  complex  frequency  response  at  spatial  location 
r  of  receive  transducer.   Also  referred  to  as 
the  complex  receive  aperture. 

DR(n*6.)  -  receive  far-field  directivity  function  or  beam 
pattern. 

y(t,r)  -  output  electrical  signal  from  receive  electro- 
acoustic  transducer  at  time  t  and  spatial  loca- 
tion r  of  transducer. 

Y(n,x)  "  frequency  (n)  and  angular  iy)    spectrum  of  output 
electrical  signal. 

As  was  mentioned  previously,  we  can  represent  the  ocean 
medium  as  a  linear,  time-variant,  space-variant,  random  filter 
The  term  "time-variant"  implies  motion  amongst  targets,  the 
ocean  surface,  discrete  point  scatterers,  and  the  transmit  and 
receive  apertures  (arrays) .   Discrete  point  scatterers  in  the 


ocean  may  include,  for  example,  gas  bubbles,  fish,  and  other 
particulate  matter.   The  time-variant  property  results  in 
both  Doppler  spread  and  spread  in  round-trip  time  delay  values. 
If  the  filter  is  time-invariant,  then  no  motion  is  implied. 
As  a  result,  there  will  be  no  Doppler  spread  and  no  spread  in 
round-trip  time  delay. 

The  term  "space-variant"  implies  that  the  sound  speed 
profile  (index  of  refraction)  of  the  ocean  is  a  function  of 
position.   The  space-variant  property  results  in  scatter  or 
angular  spread  due  to  refraction.   If  the  filter  is  space- 
invariant,  then  an  isospeed  medium  is  implied.   As  a  result, 
there  will  be  no  refraction,  and  hence,  no  scatter  or  angular 
spread  since  the  sound  rays  will  be  travelling  in  straight 
lines . 

In  addition,  since  any  motion  and/or  the  index  of  refrac- 
tion can  be  decomposed  into  a  sum  of  deterministic  (average) 
and  random  (fluctuating)  components,  these  random  components 
can  be  accounted  for  via  a  random  filter  representation  vis-a- 
vis a  deterministic  filter  representation. 

By  using  a  systems  theory  approach,  surface,  volume,  and/or 
bottom  reverberation  returns  can  be  modelled  as  the  outputs 
from  linear  filters.   In  addition,  target  returns  can  also  be 
modelled  as  filter  outputs.   Furthermore,  different  transmit 
signals  and  transmit  and  receive  directivity  functions  can 
easily  be  coupled  to  various  models  (i.e.,  transfer  functions) 
of  the  random,  inhomogeneous  ocean  medium  in  a  straightforward 


and  logical  fashion  in  order  to  study  their  effects  on  target 
detection  or  parameter  estimation  using  various  space-time 
signal  processing  algorithms. 

The  approach  of  treating  the  ocean  as  an  isospeed  medium, 
and  hence,  as  a  linear,  time-variant,  random  communication 
channel  is  well  established  [1-19].   This  linear,  time-varying, 
random  system  theory  approach  has  also  been  applied  to  target 
scattering  problems  in  radar  astronomy  [20]  and  to  communica- 
tion channels  in  general  [21-23].   However,  with  respect  to 
target  models,  past  research  efforts  have  been  devoted  mainly 
to  the  slowly  fluctuating  point  target  problem  [24-31] .   Ef- 
forts to  treat  more  complicated  target  models  were  made  by 
Kooij  [32],  Moose  [10],  and  Ziomek  and  Sibul  [19,33].   Kooij 
[32]  and  Moose  [10]  both  modelled  the  target  as  a  linear,  time- 
invariant,  deterministic  filter  while  Ziomek  and  Sibul  [19,33] 
modelled  the  target  as  a  linear,  time-varying ,  random  filter. 
In  addition,  Ziomek  [34]  has  shown  that  the  form  of  the  gen- 
eralized ambiguity  function  can  be  derived  by  treating  the 
scattered  acoustic  pressure  field  from  a  point  target  (in 
relative  motion  with  respect  to  a  bistatic  transmit/receive 
array  geometry)  as  the  output  of  a  linear,  time-varying,  random 
filter. 

Some  work  has  been  done  in  treating  the  ocean  medium  as  a 
linear,  time-variant,  space-variant,  random  filter  by  Laval 
[9,35]  and  Laval  and  Labasque  [36].   However,  the  notation 
used  to  incorporate  the  space-variant  property  is  ad  hoc,  i.e., 
spatial  variables  are  simply  included  in  the  arguments  of  the 


impulse  response  and  transfer  functions,  for  example,  rather 
than  having  eyolved  from  a  systematic  and  consistent  notation 
based  upon  linear,  time-varying,  space-varying  system  theory. 
In  addition,  Laval  and  Labasque  [36]  assume  functional  forms 
for  the  ocean  transfer  function  instead  of  deriving  them. 
Middleton  [37,38]  also  studied  underwater  acoustic  propagation 
in  a  random,  inhomogeneous  ocean,  but  did  not  concern  himself 
directly  with  the  derivation  of  random,  time-variant,  space- 
variant  ocean  transfer  functions.  He  described  the  propagation 
phenomena  using  space-time  operators. 

In  this  paper  we  will  study  underwater  acoustic  propagation 
in  a  random,  inhomogeneous  ocean  by  treating  the  ocean  medium 
as  a  linear,  time-variant,  space-variant,  random  filter.  Sec- 
tion II  is  devoted  to  a  discussion  of  the  fundamentals  of 
linear,  time-variant,  space-variant  filters  and  is  based  upon 
the  generalization  of  the  results  contained  in  Ziomek  [39] . 
A  consistent  notation  is  developed  in  a  systematic  manner.   Va- 
rious system  functions  are  introduced  and  important  input-output 
relations  and  multidimensional  (time-space)  Fourier  transform 
pairs  are  derived  for  both  deterministic  and  random  filters. 
To  the  best  of  the  author's  knowledge,  the  expressions  presented 
in  Section  II  have  not  appeared  previously  in  the  literature. 

The  equations  necessary  to  couple  the  medium's  transfer 
function  to  the  far-field  beam  patterns  of  the  transmit  and 
receive  apertures  (arrays)  and  to  the  frequency  spectrum  of 
the  transmitted  signal  are  discussed  in  Section  III. 


Section  4.1  is  devoted  to  the  main  problem  of  the  paper 
which  is  the  derivation  of  a  random  ocean  transfer  function 
incorporating  the  W.K.B.  approximation.   In  the  process  of 
the  derivation,  the  index  of  refraction  is  decomposed  into 
deterministic  and  random  components  and  it  is  assumed  that 
the  medium  is  "weakly  scattering" .   The  transfer  function 
derivation  was  motivated  by  the  work  of  Clarke  [40]  on  the 
application  of  the  W.K.B.  approximation. 

Section  4.2  is  devoted  to  an  example  calculation  of  the 
equations  for  the  random,  output  electrical  signals  appearing 
at  each  element  in  a  receive  planar  array  of  complex  weighted 
point  sources  in  terms  of  the  frequency  spectrum  of  the  trans- 
mitted electrical  signal,  the  transmit  and  receive  arrays, and 
the  random  ocean  medium  transfer  function. 

Finally,  in  Section  4.3,  the  autocorrelation  function  of 
the  transfer  function,  which  is  also  known  as  the  coherence 
function,  is  calculated. 
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II.     FUNDAMENTALS  OF  LINEAR  TIME-VARIANT  SPACE-VARIANT  FILTERS 
2. 1   Deterministic  Filters 

2.1.1   Impulse  response  and  transfer  functions 
A  linear,  time-variant,  space-variant,  filter  is 
depicted  in  Fig.  3,      where  it  is  characterized  by  its  corres- 
ponding time-varying,  space-varying  impulse  response  h(x,r  ; 
t,r)  .   The  function  h(x,r  ;  t,r)  describes  the  response  of  the 
filter  at  time  t  and  spatial  location  r  =  (x,y,z)  due  to  the 
application  of  an  unit  impulse  at  time  (t-i) ,  or  t  seconds  ago, 

and  at  a  distance  |r-r   meters  away  where  r   =  (x  ,y  ,z  ) . 

'—  — o '  *  — o     o   o   o 

Note  that 


h(x,r  ;  t,r)  =  h(t,r;  t-r,  r-r  ) .      (2.1-1) 
— o    —        —       —  — o 


The  relationship  between  the  input  signal  x(t,r) 
and  the  output  signal  y(t,r)  is  given  by 


00     oo 


-//■ 


y(t,r)  =   /   /  x(t-x,  r-r  )h(r,r  ;  t,r)dxdr 

—  — o     — o    —    — o 


-co   — »oo 


(2.1-2) 

where  it  should  be  noted  that  both  the  input  and  output  signals 
are  functions  of  time  and  space. 

Example  2.1-1 

■*- 

Different  forms  of  Eq.  (2.1-2)  can  be  obtained  by  making 
the  following  simplifying  assumptions: 
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x(  t,  x  ) 


►  y  (  t,  r  ) 


Fig.  3.   Representation  of  a  linear, 
time-variant,  space-variant, 
filter. 
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(1)   if  the  linear  filter  h  is  time- invariant  and  space-inva- 
riant, then 


h(T/ro;  t,r)  =  h(t  -  [t-c]  ,  r  -  [r-r^] ) 


=  hd,^)  (2.1-3) 


and  as  a  result,  Eq.  (2.1-2)  reduces  to 


OO     00 


y(t,r)  =   /   |x(t-T,  r-r  )h(x,r  )dtdr 

—  — o     — o    — o 


—  oo  —oo 


(2.1-4) 

which  is  a  multidimensional  convolution  integral  as  would 
be  expected. 
(2)   if  h  is  time-invariant  but  space-variant,  then 


htt,^;  t,r)  =  h(t  -  [t-x],  r;  r-r^) 


=  h ( t , r ;  £-£q  ) 


=  h(x,r  ;  r)  (2.1-5) 

— o  — 


and  as  a  result,  Eq.  (2.1-2)  reduces  to 


oo     oo 


y(t,r)  =   /   |x(t-T,  r-r  )h(x,r  ;  r)dxdr 
—      /   /        —  — o     — o  —    — o 


—  oo   — oo 


(2.1-6) 
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(3)   if  h  is  time-variant  but  space-invariant,  then 


h(r,r:  t,r)  =  h(t,  r  -  [r-r  ] ;  t-x) 
— o    —         —    —  — o 


h(t,r  ;  t~c) 
— o 


h(r/ro;  t)  (2.1-7) 


and  as  a  result,  Eq.  (2.1-2)  reduces  to 


OS     » 

r)  =   /   /  x(t-t ,  r-3 

_      j  j 


y(t,r)  =   /   /  x(t-x,  r-r^JMx,^;  tjdxdr^. 


—  oo   — oo 


(2.1-8) 
***** 

Note   that   if   Eq.     (2.1-2)    is   rewritten  as 


oo  as 

y(t,r)    =       /      Jx(t-a,    r-c.)h(a,_£;    t,r)dad£  (2.1-9) 


//■ 


-co      —  ao 


and  if  the  input  is  an  unit  impulse  applied  at  time  a  =  (t-x) 
seconds  at  a  distance  |cj  =  | r-r  j  meters  away,  i.e.,  if 

x(a,£)  =  5(a  -  [t-x],  z    -  [r-J^] ) ,     (2.1-10) 

then  substituting  Eq.  (2.1-10)  into  Eq.  (2.1-9)  yields 
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OO      GO 


*.// 


y(t,r)  =   /   /  fiCx-a,  ^"OhCa/^;  t,r)dctd£ 


—  OO   —  oo> 


=  h(x,r  ;  t,r) .  (2.1-11) 

— o    — 


Analogous  to  the  frequency  response  or  transfer 
function  H(f)  of  linear,  time- invariant  systems  is  the  time- 
varying,  space-varying,  frequency  response  or  transfer  function 
H(f,v_;  t,r)  of  linear,  time-variant,  space-variant,  systems. 
It  is  defined  as  follows: 


H(f,v;  t,r)  =  F  F   {h(x,r;  t,r) }        (2.1-12) 
—    —     x  r_     — o    — 


or 


OO      OO 


'-' '-  /  /* 


H(f,v_;    t,r)    =11  h(x,r    ;    t,  r)  exp  (-  j  2TTf  x  )  exp  (  +  j  27Tv_«r    )  dxdr 


—  OO    —  -D 


(2.1-13) 

where  f  corresponds  to  input  frequencies  in  HZ.  and  v_  is  a 
vector  whose  components  are  input  spatial  frequencies  with  units 
of  cycles/meter.   As  was  previously  mentioned  in  Section  1, 
spatial  frequencies  are  related  to  wavenumber  components,  and 
hence,  they  represent  directions  of  wave  propagation.  Similarly, 
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h(x,r    ;    t,r)    =   F~1f"1{  HCf  ,v  ;    t,r)}  (2.1-14) 

— o  —  z     v  —  — 


or 


OO  OO 


r    ;    t,r)    =       /       /  H(f,v«    t,r)  exp  (  +  j2tt  ft  )  exp  (-  j27rv_*  r    )dfdv_. 


h(T  , 

— V  —  fl  —  —  — 

-co       — ao 

(2.1-15) 

The  choice  of  a  plus  (+)  sign  in  the  exponent  of 

exp  (  +  j2ttv_«  r)  appearing  in  Eq.  (2.1-13),  which  corresponds  to 

the  forward  spatial  Fourier  transform  w.r.t.   r  ,  was  not 

— o      

arbitrary.   This  choice  of  sign  convention  is  meant  to  be  con- 
sistent with  that  of  the  spatial  Fourier  transform  relation- 
ship between  a  complex  aperture  function  and  its  directivity 
function  (beam  pattern)  as  is  developed  later  in  Section  3. 
Besides,  the  integrand  term 


exp  (- j2iTf  x  )  exp  (  +  j2irv_«  r  ) 


appearing  in  Eq .  (2.1-13)  has  the  nice  physical  interpretation 
of  being  a  time-harmonic  plane  wave  travelling  in  the  direction 
of  increasing  r   =  |r  |  with  the  sign  convention  as  given. 

Example  2.1-2 

If  the  linear  filter  h  is  time- invariant  and  space-in- 
variant, then  its  corresponding  transfer  function  can  be  ob- 
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tained  by  substituting  Eq.  C2.1-3)  into  Eq.  C2.1-13).   Doing 
so  yields  the  following  result; 


00      00 


//• 


H(f,\0  =   /   /  hCx,r  )exp(-j2TTfx)exp(+j2Trv_»ro)dxdro  . 


—  OO    — QO 


(2.1-16) 


Similarly,  from  Eq.  (2.1-15)/ 


ao  oo 


■//' 

_  oo         — OO 


h(x,ro)    =        I       I  H  (f  f\j)  exp  (  +  j2if  x)  exp  (-j2irv_«r    )dfdv_. 

(2.1-17) 


***** 


Calculating  the  transfer  function  according  to 
Eq.  (2.1-13)  assumes  that  one  knows  the  impulse  response  func- 
tion.  However,  even  if  the  impulse  response  is  not  known, 
the  transfer  function  can  still  be  obtained  by  using  a  time- 
harmonic  plane  wave  as  an  input  signal.   This  can  be  shown  by 
representing  a  linear,  time-varying,  space-varying, filter  by 
the  linear  operator  L(»)  which  operates  on  input  signals  that 
are  functions  of  both  time  and  space.   The  output  of  the 
filter  y(t,r)  can  then  be  expressed  as 

y(t,r)  =  L[x(t,r)J.  (2.1-18) 
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If  we  let 

x(t,r)  =  exp(+j2TTft)exp(-j2™-r)  ,   (2.1-19) 

which  is  in  the  form  of  a  time- harmonic  plane  wave  travelling 
in  the  direction  of  increasing  r  =  |rj  ,  then  from  Eqs.  (2.1-2! 
and  (2.1-18)  we  can  write 


OO      00 


L  [exp(  +  j2-rrf  t)  exp  (-j2ttv_«  r)  ]  /      /  exp[  +  j 2tt f  ( t-r)]  exp [-  j 2t;v_*  (r-r    )] 


-co    —  -ju 


h(x,r_;  t,r)dxdr  , 
— o    —    — o 

(2.1-20) 

and  by  using  the  definition  of  the  transfer  function  given  by 
Eq.  (2.1-13),  Eq.  (2.1-20)  can  be  expressed  as 

L[exp  (  +  j  2^f  t)  exp  (- j  27Tv_«r  )]  =  H(f,v_;  t,  r)  exp  (  +  j2^f  t)  exp  (-  j2irv_«r) 

(2.1-21) 

where  L  [exp  (  +  j  2-rrf  t)  exp  (-j2irv_»r)  ]  is  the  response  of  the  filter 
to  exp  (  +  j  2irf  t)  exp  (-  j2tt_v  «r),  and  H(f,_u;  t,r)  is  the  time-varying 
space-varying  transfer  function  of  the  filter  evaluated  at 
f  and  _v. 

Equation  (2.1-21)  is  a  fundamental  result  that 
will  be  used  in  the  derivation  of  an  ocean  transfer  function 
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in  Section  4.1.   In  addition,  Eq.  (.2.1-21)  can  be  used  to 
represent  the  output  y(t,r)  in  terras  of  the  transfer  function. 
For  example,  if  we  express  the  input  x(t,r)  in  terms  of  the 
following  multidimensional  inverse  Fourier  transform 


x(t,r)    =   F'V^XCfjv)  }  (2.1-22) 


or 


00  OO 


x(t,r)    =        J       /  X(f ,  v_)exp(  +  j2Trft)  exp(-j2TTv_«r)dfdv_, 

— '00       — oo 

(2.1-23) 

then  from  Eqs.  (2.1-18)  and  (2.1-21),  we  can  write  the  output 
y(t,r)  as 


oo  oo 


II 

—  CO       —  oo 


y(t,r)    =        /      /  X(f ,  v_)L[exp(  +  j2iTf  t)exp(-j2TTv«r)  ]dfdv_. 

(2.1-24) 


or 


OO  CO 


y(t,r)    =       /       /  X(f ,  v_)H(f  ,w    t,r)exp(  +  j2Trft)exp(-j2TTv-r)dfd: 

—   CO  —  vJ3 

(2.1-25) 
Compare  Eq.  (2.1-25)  with  Eq.  (2.1-2). 
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2.1.2   Spreading  and  bi-frequency  functions 
Two  additional  filter  functions  will  now  be  intro- 
duced; namely,  the  spreading  function  and  the  bi-frequency 
function.   The  spreading  function  will  be  discussed  first. 

The  spreading  function   S(x,r  ;  <$>,<)    is  defined  as 


SCt,^;  <j>,<_)  -  FtFr{h(T/ro;  t,r)  }       (2.1-26) 


or 


— o 


//' 


S(x,r    ;    $,k)    =        /      /Mx,r    ;    t,  r)  exp  (-  j2Tr<j>t)  exp  (  +  j  2tt  <_♦  r)  dtdr 


—  oo      —  oo 

(2.1-27) 


where  4>  corresponds  to  the  rate  of  change  of  the  filter's  im- 
pulse response  in  HZ.  and  <_   is  a  vector  whose  components  are 
spatial  frequencies  which  correspond  to  the  rate  of  change  of 
the  impulse  response  in  cycles/m.   Similarly/ 


h(r,r    ;    t,r)    =   F~1f"1{S (t , r    ;     $,<)}  (2.1-28) 

— o  —  <p      <  — o  — 


or 


oo  oo 


Mx,r    ;    t,r)    =       /      /    S(t,r    ;    <p  ,<_)  exp  (  +  j2ir<|>t)  exp  (-  j  2tt<_«£)  d*d<_. 


_oo      — oo 


(2.1-29) 
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An  alternate  representation  of  the  output  y(t,r) 
can  be  obtained  by  substituting  Eq.  (2.1-29)  into  Eq.  (2.1-2). 
Doing  so  yields 

00  OO  00  OD 


////■•" 


y(t,r)    =1111  x(t-t,    r-r    )  exp  (  +  j  2ir<J>t)  exp  (-  j2™-r) 


—  ao      — oo      — oo      — oo 


S(x,r    ;    4>,<)d<J>dKdTdr  (2.1-30) 

— o  —         —       — o 


where  the  integrand  terra 


S(T,r  ;  *,<)x(t-T.  r-r  )  exp  (+j2ir$t)  exp  (-  j  2^<  •  r) 
— o    —        —  — o  —  — 


is  a  time,  space,  and  "frequency"  (frequency  in  HZ.  and  spatial 
frequencies  in  cycles/m. )  shifted  replica  of  the  input  signal 
x(t,r)  weighted  by  the  spreading  function   S(t,r  ;  <fr,0- 
Therefore,  for  a  given  input  (source)  location  as  specified  by 
r  ,  the  spreading  function  determines  the  amount  of  spread  in 
round-trip  time  delay  t  (sec.)  ,  frequency  <J>(HZ.),  and  spatial 
frequencies  <_  (cycles/m.)  that  an  input  signal  will  undergo  as 
it  passes  through  a  linear,  time-varying,  space-varying,  channel. 
Equations  (2.1-30)  and  (2.1-2)  both  indicate  that  the  output 
y(t,r)  at  some  receiver  location  r  is  dependent  upon  the  loca- 
tion r   of  the  input  (source)  via  the  spreading  function 
S(x,r  ;  <b  ,k)    w.r.t.  Eq.  (2.1-30),  and  via  the  impulse  response 
h(x,r  ;  t,r)  w.r.t.  Eq.  (2.1-2).   For  example,  the  existence  and 
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extent  of  SOFAR  transmission  depends  on  both  the  depths  of  the 
source  and  receiver  [41] . 

The  last  filter  function  to  be  discussed  is  the 
bi-frequency  function.   The  bi-frequency  function  B(f,v_;  <j>  ,jc) 
is  defined  as 


B(f,v;  *,<)     =  FtFr{H(f/Ni;  t,r)  }  (2.1-31) 


or 


oo    oo 


B(f,v_;  ♦/£)  =   /   1  H(f  ,v_;  t,  r)  exp  (-  j27rcj>t)  exp  (  +  j  2tt<_-  r)  dtdr . 


■oo   — oo 


(2.1-32) 
Similarly, 


H(f,v;  t.r)    =  F^1F~1{B(f,v;  +,  ic)  }  (2.1-33 


or 


OO      00 


H(f,v_;  t,r)  =   /   /  B(f,u_;  $  ,  <)  exp  (+j2ir$t)  exp  (-j2ttk_'£)  d<j>d<_. 

—  CO       —  oo 

(2.1-34) 

The  term  bi-frequency  function  was  originally  used  to  denote 
the  appearance  of  the  two  frequency  variables  f  and  <j>  in  linear, 
time-varying  systems  theory.   We  can  generalize  this  notion 
by  using  the  term  bi-frequency  function  to  denote  the  appearance 
of  the  input  "frequency"  pair  f  and  v,and  the  filter  variation 
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"frequency"  pair  <j>  and  <_   in  linear,  time-varying,  space-vary- 
ing, systems  theory. 

Just  as  the  spreading  function  gives  an  indication 
of  how  rapidly  the  impulse  response  changes  with  time  and 
space,  the  bi-frequency  function  gives  an  indication  of  how 
rapidly  the  transfer  function  changes  with  time  and  space. 

The  bi-frequency  function  can  also  be  obtained  by 

taking  the  Fourier  transform  of  the  spreading  function  w.r.t. 

t  and  r  ,  i.e., 
— o 


B(f,v;  <*>,<)  =  FTFr  {S(T,rQ;  <j>,<_)}        (2.1-35) 

— o 


or 


CO      CO 


// 


B(f,y_;    <}>,£)    =       /      /    S(x,r    ;    <£  ,  <_)  exp  (-  j  2irf  x  )  exp  (  +  j  27r\j_«r    )dxdr 


—  CO        —  CO 


(2.1-36) 


and ,    s  imi lar ly , 


Sfx,^;    <*>,<)    =    F^1F~1{B(f,v;     ?,<)}  (2.1-37) 


or 


S(x,r    ;    4>,k)    =      /       I   B(f,v;     <p ,  <  )  exp  (  +  j  2-nf  x  )  exp  (-  j  2-au  «r    )dfdv    . 

-         j     j  -        -  --a        - 

—  CO       —  oo 

(2.1-38) 
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In  addition,  the  impulse  response  and  bi-frequency 
functions  form  a  Fourier  transform  pair,  i.e., 


B(f,v;    <*>,<)    =    FTFr   FtFr{h(x,ro;    t,r)  }  (2.1-39) 


or 


CO  GO  OO  'JO 


B 


(f,v;    $,<)    =       I       I      I       I  h^T'£0;    t,r)  exp(-j27rf  x)  exp  (  +  j2Tr_y  *r 


«>  oo         —  OO         —  OO         —CO 


exp  (-j  2ir<j>t)  exp  (  +  j2Ti<_«r)  dxdr   dtdr 

(2.1-40) 
and,    similarly, 


Mt,^;    t,r)    =    F^1F^1F^1F^1{B(f,v;     f,<)} 


(2.1-41; 


or 


OO 


OO  OO  OO 


hft,^;    t,r)    =       I       I      I      I   B(f,v_;    <|»  ,<_)  exp  (  +  j  2uf  x)  exp  (-  j2ttv  •  r^) 


—  OO         —  OO         —  OO         —OO 


exp  (  +  j2ir<j>t)  exp  (-j2irjc-r)  df  dvd<j)d<_. 

(2.1-42) 
The  interdependence  which  exists  amongst  the  four 
filter  functions  is  illustrated  in  Fig.  4.       Any  one  of 
these  functions  may  be  used  to  define  completely  a  linear, 
time-varying,  space-varying,  communication  channel. 
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TIME- VARYING  SPACE- VARYING 
IMPULSE  RESPONSE 


SPREADING 
FUNCTION 


H(f,*;  t,   r  ) 


TIME-VARYING 
SPACE -VARYING 
TRANSFER 
FUNCTION 


B  Cf,  j? ;  4>,  k  ) 


BI-FREOUENCY 
FUNCTION 


Fig.  4.   Interdependence  amongst  the  four  filter 
functions  that  characterize  linear,  time- 
varying,  space-varying,  communication 
channels. 
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2.1.3   Output  frequency  and  angular  spectrum 
The  output  frequency  (n )  and  angular  C0_)  spectrum 
Y(n,S)  is  defined  as 


Y(n,S)  =  FtFrCy(t,r)}  (2.1-43) 


or 


00      00 


"■•-' »  /  / 


Y(n,B_)  =   /   /  y(t,r)exp(-j2Trnt)exp(  +  j2Tr8_«r)dtdr     (2.1-44) 


_00    _00 


where  1  corresponds  to  output  frequencies  in  HZ .  and  8  is  a 
vector  whose  components  are  output  spatial  frequencies  with 
units  of  cycles/meter.   Substituting  Eq.  (2.1-25)  into  Eq. 
(2.1-44)  yields 


»■// 


CO  00  oo  CO 


Y(n,S_)    =      /      /   X(f,v)       /      /   H(f,v;    t,r)exp[-j2*(n-f)t]' 


•  OO        M   -O  —  oo        «•  oo 

exp  [+j2ir  (8-v)  «r]dtdrdfdv 


(2.1-45) 

where,  from  Eq.  (2.1-32),  it  can  be  seen  that  the  inner  multi- 
dimensional integral  is  equal  to  B(f,_u;  n-f,  S_-v_)  so  that 

oo     oo 

Y(n,6_)  =   /   /  X(f,v)B(f,v;  n-f,  6_- v )  df  dv 

—  oo   — oo 

(2.1-46) 
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which  is  in  the  form  of  a  multidimensional   "frequency"  domain 
convolution  integral.   Note  that  <j>  =  n-f  or  n  =    f+<t>,  and  that 
<_  =    8_-v_  or  £  =  \H-£.   Thus,  the  output  frequencies  n  in  HZ.  are 
equal  to  the  sum  of  the  input  frequencies  f  and  the  variations 
in  frequency,  <J> ,  due  to  the  time-varying  property  of  the  filter, 
Similarly,  the  output  spatial  frequencies  _S  (directions  of 
wave  propagation)  in  cycles/m.  are  equal  to  the  sum  of  the  in- 
put spatial  frequencies  v_  (directions  of  propagation  of  trans- 
mitted plane  waves)  and  the  variations  in  spatial  frequency 
(directions  of  wave  propagation)  ,  <_,    due  to  the  space-varying 
property  of  the  filter.   Equation  (2.1-46)  demonstrates  that  a 
linear,  time-variant,  space-variant,  filter  will  spread  the 
input  frequency  and  angular  spectrum  X(f,v_)  in  both  frequencies 
in  HZ.  and  spatial  frequencies  in  cycles/m. 

Example  2.1-3 

If  the  linear  filter  h  is  time- invariant  and  space-in- 
variant, then  H(f,v_;  t,r)=  H(f,v_)  [see  Eq.  (2.1-16)]  and,  as  a 
result,  Eq.  (2.1-32)  reduces  to 


B(f,v;     <t> ,  <)    =   H(f,v)        /   exp(-j2Tr<$t)dt       /   exp  (  +  j  2  tt  <•  r )  dr 


or 


B(ffv;     <Ji,k)     =    H(f  ,v)  6(<|>)  5(k)  (2.1-47) 
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since 


and 


FtU>  =  5(<j>) 


Fr(l}  =  S(<) 


Substituting  Eq.  (2.1-47)  into  Eq.  (2.1-46)  yields 


00  00 


"">-  // 


Y(n,B)    =       /      /   X(f ,v)H(f , v)6 (n-f) 6 (S-v)dfdv 


—  oo       *>oo 


which   simplifies    to 

Y(n,£)    -  X(n#e.)H(n,3.) , 

and  by  replacing  n  with  f  and  _6_  with  v_,    we  finally  obtain 

Y{£,\>)    =  X(f  ,v)H(f  ,v_)  (2.1-48) 

which  is  the  output  frequency  and  angular  spectrum  from  a  linear, 
time- invariant ,  space- invariant,  filter.   Note  that  the  output 
frequency  f  and  the  output  spatial  frequency  vector  \>_   are  iden- 
tical with  those  of  the  input  X(.£,v).      Hence,  as  would  be  ex- 
pected, there  is  no  "frequency"  spreading.   This  would  corres- 
pond to  the  physical  situation  of  transmitting  a  signal  via  a 
transmit  aperture  (array)  to  a  receive  aperture  (array)  when 
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the  platforms  containing  the  apertures  (arrays)  are  not  in 
motion,  and  the  intervening  ocean  medium  has  a  constant  speed 
of  sound  (index  of  refraction)  and  no  discrete  point  scatterers 
in  motion.   No  motion  implies  no  Doppler  (frequency)  spread 
and  a  constant  speed  of  sound  implies  that  the  sound  rays 
travel  in  straight  lines,  i.e.,  there  is  no  angular  spread 
(scatter)  or  change  in  the  direction  of  propagation  of  the 
transmitted  sound  rays. 
***** 
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2 . 2   Random  Filters 

2.2.1   Filter  autocorrelation  functions 
In  Section  2.1,  four  system  functions  were  intro- 
duced which  are  used  to  characterize  linear,  time-variant, 
space-variant,  filters.   However,  if  the  filter  is  random, 
then  each  of  these  system  functions  must  be  considered  as 
random  functions.   As  a  result,  we  must  work  with  the  respec- 
tive system  autocorrelation  functions  which  are  defined  as 
follows : 


rk<t't,'£^'£JI'*  t,t\r,r')  =E{h(r,r  ;  t,r)h*  (t  » ,r! ;  t\r')} 

n      — o  - o       —  —  — o    —       — o     — 

(2.2-1) 


Rjj(f  ,f  '  ,v,v' ;  t,t',£,£*)  =  E{H(f,v;t,r)H*(f  ,,v';  t\r')} 

(2.2-2) 


RQ  (t,t  '  ,r  ,r!  ;  $,$',<,<•)  =  E{S(x,r  ;  $ ,  tc)  S*  (t  '  ,r*  ;  **,<•)} 

o         —^j     — O  —  — -  - "O      "■■  — -O 

(2.2-3) 
and 

RgCf/f  ,v,v';  ♦,♦', £.,£.')  =  E{B(f,v;  $,k)B*  (f '  ,  v  ' ;  $•,<•)} 

(2.2-4) 
where  E{ • }  is  the  expectation  or  ensemble  average  operator 
and  the  asterisk  denotes  complex  conjugate. 

If  we  use  the  sign  convention  that  forward  trans- 
forms w.r.t.  x,t,r',  and  r*  are  defined  with  a  minus  sign  in 
the  exponent  of  the  complex  exponential  ( inverse  transforms 
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w.r.t.  f,<f>,v_' ,  and  k_%    are  defined  with  a  plus  sign)  and  for- 
ward transforms  w.r.t.  t1, t',r  ,  and  r  are  defined  with  a 
plus  sign  (inverse  transforms  w.r.t.  f',^,^,  and  <_   are  defined 
with  a  minus  sign) ,  then  it  can  be  shown  that  the  four  system 
autocorrelation  functions  are  related  by  the  following  Fourier 
transform  pairs: 


R^fff  \v,v\-  t,t\  £,£')  =  F  FT,Fr  Fr ,  {1^  (t  ,  x  '  ,rQ,  r^;  t,t',r,r' 

— o  — o 

(2.2-5) 

and 

Rhl*'*,'Z0'*i''    t,t\r,r')  =  P^1P^P^1P^{RH(f,f*fvfv,s  t.t'.r,!1 

(2.2-6) 


R^fff  ',v,v';  <J>,4>  '/</<')  =FtFt-FrFr,{RHtf'f,'-'-,;  t,f ,£,£*)} 

(2.2-7) 

and 

RjjCf/f*  rVi/V'jt^f  ,£,£')  =  F~1F^F^1F^{RB(f  ,f  ' , v, v  '  ;  *  ,*',<,<_')}  ; 

(2.2-8) 

rS(t't,'£o'^o;    *'*''J£'£,)    =  FtFt,FrFr'  {Rh(T'T'  'r^'r^;    fc't,'£'£,)| 

(2.2-9) 

and 

R.   (t,t'  ,r    ,rl;     t,f  ,r,r')    =    f'V^F      f'.CRq  (t  ,  t  r,r    ,£•  ;     $  ,<J>',k,k  ' 
a  — O    — O  —   —  <J>        <j>        K       J£  ^  — O    — O  —   — 

(2. 2-10) 


R^fjf'^v';     4>    *',£,£*)    =   F   F      F      P.fHJT.r'.r    ,£■;     <f> ,  *  '  ,  k  ,  k  •  )'* 

— o  — o 

(2.2-11) 
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and 

(2.2-12) 
and  finally, 

Rgff^'^v''  Kf.M')  =  FTFT,Fr  Fr,FtFt,FrFr,{Rh(T/T'/ro/r<|); 

— o  — o      —  — 
t,fr,r')}  (2.2-13) 

and 


VT'T''£o'^>;  t't\£fr')  =  F-1F-}F-1F-^1F;}F;1F2^{RB(f,f  ,v,v', 

*,<j>' ,<,<*)  K  (2.2-14) 


The  interdependence  which  exists  amongst  the  four 
filter  autocorrelation  functions  is  illustrated  in  Fig.  5. 
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2.2.2   Uncorrelated  spreading  and  the  scattering 
function 

In  this  section  we  will  examine  the  consequences 
of  assuming  that  the  spreading  function  S(x,r  ;  <$>,<)    is  un- 
correlated with  S(t',  r '  ;  i>'/£')  for  all  values  of  t'  ^    t, 

r  *  ?*  r  ,  & '  ?    <p ,    and  < '  ^  < . 
— o   — o  —    — 

The  assumption  of  uncorrelated  spreading  is  mathe- 
matically equivalent  to  stating  that  the  autocovariance  of 

S(x,r  ;<J>,kc)  and  S(T',r';  <j> '  ,  < ' )  is  zero  for  all  values  of 
— o    —  — o      —      — — -     — — " 

t*  ?    t,  r^^r^,  <j>'  ?    0,  and  <_%    ±    <_,    i.e., 


cS(t'  t'  '£o'^o;  *'  *'  '-'-'*  =  Rs(T/T*  ,ro'^o;  *'*''£'£')  ~ 

t'  ^  t,  r^^r^,  <f>'  f    cp,  and  £*  ^  £ 

(2.2-15) 


where  Cc  is  the  autocovariance  function  and  yc(r,r  ;  $,<)  = 
E{S(r,r  ;  $ ,  k)   }.   If  it  is  assumed  that  yc(T,r  ;  4> ,  < )  =  0, 
then  Eq.  (2.2-15)  is  equivalent  to 


Rs  (x,  t  '  /r^^r^;  <J>,  <j>'  ,£,£'  )  =  Rg  (  t  ,  r^;  *  ,jc)  5  (  t-t  *  )  • 

5  (r  -r')  6  U-cf)')  5  (<-<')    (2.2-16) 
— o  — o  —  — 


where 


R3(T,ro;j),£)  =E{|S(x,ro;  <|. ,  <)  |  2  }  (2.2-17) 
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is  called  the  scattering  function  and  is  equal  to  the  mean 
squared  value  of  the  spreading  function.   The  scattering  func- 
tion can  be  thought  of  as  an  average  power  density  function 
which  determines  the  average  amount  of  spread  that  an  input 
signal's  power  will  undergo  as  a  function  of  round-trip  time 
delay  x  ,    frequency  $  ,  and  spatial  frequencies  <_  for   a  given 
input  (source)  location  r  .   Note  that  the  scattering  function 
is  a  real,  positive  valued  function. 

Equation  (2.2-16)  is  the  result  of  the  assumption 
that  the  spreading  function  is  zero  mean.   However,  if  the 
spreading  function  is  non-zero  mean,  it  is  convenient  to  do 
the  analysis  with  the  centered  process 


S^(x,r  ;  $,<)  =  S(i,r  ;  4> ,  < )  -  Uc(T,r  ;  *,<). 

C    — O      —  — O      —       b     — O      — 


The  random  process  S  (x,r  ;  4>,<_)  has  zero  mean,  and  as  a  re- 
sult, its  autocovariance  function  is  equal  to 


cq  'T'T''£^'£^'  4>,*,,k,i')  =  E{Sr(x,r  ;  ♦,<)  S*  (t  •  ,r» ;  *',<_')> 

=  Rs   (T,T»  /^r^  ;   $,<j>'  ,£,<') 

=  CgC?,?',^,^;  4>,4>,,<,<I)  , 


i.e.,  the  autocorrelation  function  of  the  zero  mean  centered 
process  is  equal  to  the  autocovariance  function  of  the  original 
non-zero  mean  spreading  function. 
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Let  us  next  examine  the  effect  that  the  assumption 
of  uncorrelated  spreading  has  on  the  remaining  three  system 
autocorrelation  functions.   As  Fig.  5     indicates,  R^  can 
be  obtained  from  Rq  by  performing  the  following  transformations 


Rjjtf.f.v.v-  ,t.f.r,r->.   F^F^VIf^  ,Fr   F    .    {Rgtr.f, 

—       —  — o     — O 

r^rl;    4>,<|>  '  ,<_,<_*)}  (2.2-18) 

or 


^(f  ,f 'f^,v'  ;   t,t',r,r*)   =     /     ...      /  Rs(t,t'  iXqiZ^''    <b,  <t> ',<_,<_' ) 


I ...  /, 


exp[+j2.TT  ($t-$ '  t*)]exp[-j2ir  (<_•£-<_'  •£')]• 

exp[-j2TT  (f  t-f  'x  ' )  ]exp[  +  j2TT  (v_»r   -v/  -r!)  ] 

d^dct'dKdK'dTdT'dr   dr '     .  (2.2-19) 

—  —  — o   — o 

If  Eq.  (2.2-16)  is  substituted  into  Eq.  (2.2-19),  then 


R^(f ,f ' rv»v* ;    t,t',r,r')    =RH(Af,Aw    At,Ar)       (2.2- 


20) 


where 


00  00  CO  CO 


RjjtAf^Av;    At,Ar)    =      /       /       /       /     Rgd,^;    <j> ,  k)  exp  (  +  j  2tt  <d  a  t) 


-co       — ■  co       —  aa       —  oo 


exp  ( -  j  2tt  <_•  Ar)  exp  (- j  2^Af x ) 

exp  (  +  j2ttAv  *r    )d<J>di<dxdr    , 
r      J        —  — o  —        — o 


(2.2-21) 
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Af  =  f-f'f  Av  =  v-v1 ,  At  =  t-t'f  and  Ar  =  r-r ' . 

It  can  be  seen  from  Eq.  (2.2-20)  that  when  uncorrelated  spread- 
ing is  assumed,  the  autocorrelation  function  Ru  becomes  a  func- 

rl 

tion  of  the  differences  Af,  Av,  At,  and  Ar  only.   This  implies 
that  the  random  process  H(f,v_;  t,r)  is  wide-sense  stationary 
in  frequency,  spatial  frequencies,  time,  and  space.   An  addi- 
tional requirement  for  H(f,_u;  t,r)  to  be  wide-sense  stationary 
is  that  the  mean  value 


MH(f,v;  t,r)  =  E{H(f,v;  t,r)  } 


is  a  constant.   Since  the  four  filter  functions  are  related 
by  linear  transformations  (see  Fig.  4)  ,     and  since  it  was 
assumed  that   iic(T,r  ;  $,<)    =    0,  then  M,(x,r;  t,r)  =  0, 
MH(f,v_;  t,r)  =  0,  and  nB(f,v_;  4>,<_)  =0  which  are  constant, 
zero  mean  values.   Therefore,  the  condition  of  uncorrelated 
spreading  in  round-trip  time  delay  x,  input  (source)  location 
r  ,  frequency  $ ,  and  spatial  frequencies  <_  is  equivalent  to  a 
condition  of  wide-sense  stationarity  in  frequency  Af,  spatial 
frequencies  Av_,  time  At,  and  output  (receiver)  location  Ar, 
respectively.   If  uncorrelated  spreading  in  x,  r  ,  4>,  and  <_ 
occur  together,  then  we  have  a  wide-sense  stationary  uncorre- 
lated spreading  (WSSUS)  communication  channel. 

Consider  the  autocorrelation  function  R.  next. 

n 

If  Eq.  (2.2-16)  is  substituted  into  Eq.  (2.2-10)   and  the  in- 
dicated transformations  are  performed,  then 
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Ru(T,T*/r    ,r*s    t,t*,r,r')    =R^(x/r    ;    At,  Ar)  5  (t-t  ' )  5  (r  -r » ) 
n  - o   — o  —  —  n        — o  —  — o  — o 

(2.2-22) 

where 


00  oo 


—  OS       — oo 


^(t,^;    At,Ar)    =     /       /     Rgtx,!^;    <j> ,  <_)  exp  (  +  j  2774,  A  t)  • 

exp  (-  j2ttic_»  Ar)  d<£d<_.  (2.2-23) 


Equation  (2.2-22)  indicates  that  the  random  process  h(x,r  ;  t,r) 

is  wide-sense  stationary  in  time  and  space  because  of  the  At 

and  Ar  dependence  and  since  yh(r,r  ;  t,r)  =  0.   Equation  (2.2-22) 

also  indicates  that  h(x,r  ;  t,r)  is  uncorrelated  for  all  values 
of  x '  7*  x  and  r"  7*  r  . 

Finally,  if  Eq.  (2.2-16)  is  substituted  into 
Eq.  (2.2-11)  and  the  indicated  transformations  are  performed, 
then 


Rg(f  ,f  '  ,  v,  v_'  ;    <J>,  <j> ',<_,<_' )    =   Rg(Af ,  Aw  <j> ,  <_)  6  dj)-^*  )  5  U-<_* ) 

(2.2-24) 
where 


go  co 


RB(Af,Aw     $,jc)    =      /      I     Rs(T,r    ;    $  ,  <_)  exp  (-  j  2tt  Af  x)  • 

—  00      —  oo 

exp(  +  j27TAvr  )dxdr  .  (2.2-25) 

c  —   — o    — o 

Equation  (2.2-24)  indicates  that  the  random  process  B(f,_v;  <j>,<_) 
is  wide-sense  stationary  in  frequency  and  spatial  frequencies 
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because  of  the  Af  and  Av_  dependence  and  since  y_(.f,;v;  4>,k_)  =  0. 
Equation  (2.2-24)  also  indicates  that  B(f,_v;  $,<_)  is  uncorre- 
lated  for  all  values  of  <fr '  ^  <J>  and  jc*  ^   ic. 

Therefore,  in  summary,  under  the  assumption  of 
uncorrelated  spreading  the  four  filter  autocorrelation  func- 
tions reduce  as  follows : 


6  (4>-*' )  <5  («-<■ )  (2.2-16) 


Rgtfjf  '/V,^  jt^Sr.r')    -  RjjUfjAv;    At,Ar)  (2.2-20) 

^(t,T',r    ,r!;    t,t,,r,rl)    -R^Ct^r    ;    At,  Ar)  5  (t-t  ' )  5  (r   -r« ) 
n  — o  — o  —  —  n       — o  —  — o  — o 

(2.2-22) 
and 

Rgtf  ,f  ,,v,v«;     4>,<fr»,  £,<_')     =    R^A^Av^jcMU-VJSU-^') 

(2.2-24) 


where  RQ(T,r  ;  <}>,<)    is  the  scattering  function. 
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and 


2.2.3   The  scattering  function  and  its  Fourier 

transforms 
By  inspecting  Eq.  (2.2-23),  we  can  write  that 

Rv,(T'^  '"  At,Ar)  =  F*  F~  {Rc(T,r^  ;  $,<)} 
h   — o       —     <|>  <_       s   — o     — 

(2.2-26) 

Rs(t/£o  ;  *,<)  =  *At*A£{VT'£o  '*  At'AI)} 

(2.2-27) 


or 


00     oo 


Rs  (  t  ,  r^    ;  $ ,  k)    = 


J   /  &  (T,r^  ;  At,  Ar)  exp(-j27Tcj)At)  • 


—  oo   — oo 


that 


and 


exp(+j2™«  Ar)dAtdAr.  (2.2-28) 

Also,  by  inspecting  Eq.  (2.2-25),  we  can  write 


Rq  (Af ,  Av    ;     $,*.)     =    F    Fr.  {Rs  (t,^    ;     cj> ,  <_)  } 

— o 

(2.2-29) 


Rs(x,£o    ;     4>,<)     =    F^F^(RB(Af,Av    •     ♦,<)} 

(2.2-30) 


or 


oo  oa 


RS(T'=o    ;     *'-*    =      I       I     RB(Af'A^    ;     *fl)exp(  +  j2^Afx) 


—  oo       —  oo 


exp  (-j27TA_u«r    )dAfdAv_.  (2.2-31) 


40 


Additional  transform  pairs  can  be  obtained  as 
follows.   With  the  use  of  Eq.  (2.2-23),  Eq.  (2.2-21)  can  be 
expressed  as 


RjjUfjAv    ;    At,Ar)    =   F   F      {R^x,^    ;     At,Ar) } 

— o 

(2.2-32) 

or 


oo  oo 


RjjUffAx^    ;    At,Ar)    =      J       /     R^x,^    ;    At,  Ar)  exp  (-  j2irAf  t  )• 


—  oo      — oo 


and 


exp(  +  j2TTAvr    )dxdr  (2.2-33) 

c      j        —  — o  — o 


Rh(x'^0    '"     At'A£>     =   FAfFAi{RH(Af ,A-    '"     At'A£)} 


(2.2-34) 
or 


00  oo 


// 

—  oo      — oo 


R.(x,r       ;    At,    Ar)    =11     R^(Af,Av    ;    At ,  Ar)  exp  (+j27rAf  x  )• 


exp(-j2TTAv_«r    )dAfdAv_.  (2.2-35) 


Next,    if   Eq.     (2.2-28)    is    substituted    into   Eq.     (2.2-25),  then 


RgUf/Av.    ;     <»,<)    =   FTFr   FAtFAr^Rh(x'£0    '"     At,Ar)} 

— o  — 

(2.2-36) 
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or 


CO  QO  OQ  CO 


//// 

—  OO       —CO       — oo        —  go 


^(Af^Av    ;    *,0    =      J        I       J       J     ^(t,^    ;    At ,  Ar)  exp  (-  j  2irAf  x ) 


exp  (  +  j27rA\^*  r    )  exp  (- j  2tt cj> A t) 
exp  (  +  j2n<_«  Ar)  dtdrdAtdAr 


(2.2-37) 

and 

R.   (T/r       ;    At,Ar)    =    f"Jf"1f"1f"1{ R_ (Af , Av    ;    $  ,  k  )  } 
n         — O  —  At    Av     $       <  o  —  — 

(2.2-38) 

or 

CO  CO  CO  CO 

R^d,^    ;    At,Ar)    =      J       j       I        /^(Af^Av    ;    <j> ,  <_)  exp  (  +  j  2tt  Af  t  )  • 


o  OO         •"  CO         —  CO         "■  CO 


exp  (-  j2irAv_«  r    )  exp  (  +  j 2tt<j> At)  ■ 
exp  (- j  2tt<_«  Ar )  dAf  dA\jd<pd<_. 

(2.2-39) 

Finally,  with  the  use  of  Eq.  (2.2-33),  Eq.  (2.2-37)  can  be 
expressed  as 

Rg(Af,Av  ;  <}>,<)  =  FAtFAr{RR(Af ,  A  v  ;  At,Ar)} 

(2.2-40) 
or 

»      CO 

RB(Af,Av_  ;  $,k)    =      I       I  R  (Af,A\^  ;  At ,  Ar)  exp  (-  j  2  tt  0  At)  • 


—  CO    —  QO 


exp(  +  j27nc- Ar)dAtdAr  (2.2-41) 
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and 

R^fAf.Au     :     Ah.Ar)     =    F 


RgCAf,Av    ;    At,Ar)    =   F~  F~    {R^A^A^    ;    <}>,<_)} 


(2.2-42) 
or 


QO  00 


RjjCAf/Ax^    ;    At,Ar)    =      /      /     RgCAf^Ax^    ;    <j> ,  <_)  exp  (  +  j  2Tr<f>At)  • 

—  a       — oo 

exp  (-j2tt<_»  at)  d<))d<_.  (2.2-43) 

The  various  Fourier  transform  pairs  are  summarized 
in  Fig.  6.  The  scattering  function,  or  any  of  its  Fourier 
transforms,  is  a  complete  second  order  statistical  description 
of  a  WSSUS  communication  channel. 

Another  very  important  Fourier  transform  relation- 
ship can  be  obtained  by  substituting  Eq.  (2.2-35)  into  Eq. 
(2.2-28)  which  yields 

VT'£o  ;  *'£)     =  FAtFArFAfFA^VAf'A^  ;  AtfA£>} 

(2.2-44) 

or 


00  oo  oo  oo 


S  "'^o    '     y'-'  /  H 


//// 


RJx,r       ;     $,k)    =      /       /       /       j     R    (Af,Av    ;    At,  Ar)  exp  (-j  2TT$At) 


—  oo      — oo      -»oo      —  oo 


exp  ( +  j  2  7T  <_ •  Ar )  exp  ( +  j  2 1 A f  x  )  • 
exp  (-  j2TrAv_«r  )  dAtdArdAf  dAv_ 

(2.2-45) 
and 
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RjjUfjAv    ;     At,Ar)     = 


F^F^F   F       {R-(.T,r       ;     $,<)  } 


(2.2-46) 


since 


RH(Af,Av     ;     At,Ar) 


CO  GO  CO  CO 

//// 

—  go      «■  co      «  oo      —  oo 


RgCt,^    ;    (j),  <_)  exp(  +  j2^*At) 


exp  (- j2tt<_*  Ar)  exp  (- j  2irAf  x)  • 
exp(  +  j2irAvr    )d4>d<d~dr    . 

r »     J  _   __o'  _  —o 

(2.2-21) 
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2.2.4   Input-output  relations 

We  will  now  proceed  to  derive  an  expression  for 
the  autocorrelation  function  of  the  output  from  a  linear, 
time- varying,  space- varying,  random  filter.   The  output  auto- 
correlation function  R  (t,t,,r,r' )  is  defined  as  follows: 


R  (t,f ,£,£/)  =  E{y(t',r)y*(t  r '  )  }.       (2.2-47) 


If  Eq.  (2.1-30)  is  substituted  into  Eq.  (2.2-47),  then  the 
output  correlation  function  in  time  and  space  is  given  by 


OD  00 


Ry<t't,'£'£,>  =  J  ••t(  x(t-T,  r-ro)x*(t,-Tl ,  r'-r^)- 

exp[+j2ir  Ut-<J>  *  t*)]exp[-j2ii  U-r    -   £*•£')] 

Re  (t,  t  '  ,r    ,r!    ;    $ ,  $  * ,  k,  k  *  Jd^dicdxdr    • 

d^'dic'dT'dr'     .  (2.2-48) 

—  — o 


If  the  random  filter  is  a  WSSUS  (wide-sense  stationary  uncor- 
rected spreading)  communication  channel,  then  the  autocorre- 
lation function  of  the  spreading  function  is  given  by  Eq. 
(2.2-16).   Substituting  Eq.  (2.2-16)  into  Eq.  (2.2-48)  yields 
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00  OQl  oo  to 


//// 

«GO        —  OO         —  OD>       r-  OO 


R    lt,t'  ,r,r')    =     /      /      I       /     xCt-T,r-ro)x*(t'-T,r,-ro) 


exp  (  + j  2 it <j>  At)  exp  (-  j2tt£«  Ar) 


R_(t,r      ;     (J>,  <)  d^dicdrdr    . 

b  — O  —  —  — O 


(2.2-49) 


The  mean  squared  value  of  y(t,r)  is  given  by 


E{ |y(t,r) |2}  =  R  (t,t,r,r) 


where,  from  Eq.  (2.2-48),  we  have  in  general  that 


OO  JO 

,,J  -  f-.j 


R    (t,t,r,r)    =     I   •••    /    x(t-t,r-r    ) x* (t-x ' ,r-r » ) • 
y  o  o 

exp  (  +  j2irA<J>t)  exp  (-j2tA<  «r) 


Rc(T,T',r    ,  r'    ;    $,•$* ,  k,  k'  Jd^dicdrdr    • 

d^'d^'dx'dr'  (2.2-50) 


where  A<J>  =  4>-<j>'  and  A<_  =  <_-<_'.   In  the  case  of  a  WSSUS  commu- 
nication channel,  we  have  from  Eq.  (2.2-4  9)  that 
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CO  CO  CO  CO 


yt,t,r,r)    =     f      f     J      j 


»   OO         —  OO         —  OO         —   OO 


2 
x(t-i,    r-r^)  |    R^t,^    ;     4>,£) 


dcjjd^drdr^.  (2.2-51: 


If  we  define  the  energy  of  the  output  signal  E^  as 


OO      CO 

A 


// 

—  CD   —CO 


EY  I       I     |y(t,r) |*dtdr,  (2.2-52) 


then  the  average  output  energy  E   can  be  expressed  as 


CO  CO 


_  CO        —  CO 


E      =   E{E    }    =     /      /    E{ |y(t,r) |2} 


or 

CD  CO 


E      =     /       /     R    (t,t,r,r)dtdr.  (2.2-54) 


—  ao       -co 


If   Eq.     (2.2-51)    is    substituted   into   Eq. (2.2-54),    then 


OO  OO  OO  CO 


VEx   =  R3(T'^    '     ^^d^Kdxdr^ 


—  CO  —  'JO  —JO  —CO 


(2.2-55) 
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where 


OO      CO. 


Ex  =  f  J  lxtt-£ 


I 2dtdr  (2.2-56) 


—  CO    — OO 


is  the  energy  of  the  input  signal.   Equation  (2.2-55)  indicates 
that  the  ratio  of  the  output  (received)  average  energy  to  the 
input  (transmitted)  energy  for  a  WSSUS  communication  channel 
can  be  obtained  by  integrating  the  scattering  function  of  the 
channel.   Also  note  that  the  average  output  energy  is  not  a 
function  of  the  input  signal's  shape. 

Alternate  expressions  for  the  output  autocorrela- 
tion function  can  be  obtained  from  Eq.  (2.1-25).   If  Eq. 
(2.1-25)  is  substituted  into  Eq.  (2.2-47),  then 


CO  OO  OO  OO 


t,t\r,r')    =     [     f     f     f 


R    (t,t\r,r')    =  X(ffv)X*(f,vf) 


—  OO  "» OO  ™  OO  ■■  OO 


Rjj(fff  \v,v'    ;    t,f,r,r')  • 

exp[  +  j27T  (ft-f  '  t')]exp[-j2Tr  (v-r-v_'  -r  * )  ] 
dfdvKif'dv'  .  (2.2-57) 

If  the  random  filter  is  a  WSSUS  communication  channel,  then 
the  autocorrelation  function  of  the  transfer  function  is  given 
by  Eq.  (2.2-20).   Substituting  Eq .  (2.2-20)  into  Eq.  (2.2-57) 
yields 
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CO  OO  OT  CO 


//// 


R    (t,t\r,r')    =     J      /       /       /     XCf,v)X*tf,,v,)RH(Af,Av    ;    At,Ar)- 


y 

—  GO  -»   00  —  OO  ~  0O 

exp[  +  j2Tr(ft   -    £'  t')  Jexp[-j2ir  (v»r   -v.*  •£')] 
dfdvdf'dv1.  (2.2-58) 


A  relationship  between  the  input  and  output  power  spectral 
densities  will  be  obtained  next. 

Let  us  assume  that  the  input  signal  x(t,r)  is  a 
zero-mean,  wide-sense  stationary  (in  time  and  space) ,  random 
process  which  is  uncorrelated  with  the  transfer  function. 
Under  these  assumptions,  the  output  autocorrelation  function 
given  by  Eq.  (2-2-57)  becomes 


oo    oo    oo    oo 


R  (t,f  ,r,r')  =  J   |  f     I     E{X(f,v)X*(f  ,,v')  }• 


—  OO    m»  CO    —CO   «  CO 


RH(f,f',v,v'  ;  t,f  ,£,£*)• 

exp[+j27r  (ft-f  ,t,)]exp[-j27r  (v/r-v.1  -rf )  ] 
dfdvdf'dv'  .  (2.2-59) 

Since  x(t,r)  was  assumed  to  be  wide-sense  stationary  in  both 
time  and  space,  then  it  can  be  shown  that 


E{X(f,v)X*(f,v') }  =  S  (f ,v)o(f-f ') o(v-v') 

(2.2-60) 
where 

X(f  ,v)  =  FtFr(x(t,r) }, 
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Sx(f,v)  4  FAfcFAr{Rx(At,Ar) .  }■  (2.2-61) 


is  the  power  spectral  density  of  the  input,  and 


R  (At,Ar)  =  E{x(t,r)x* (f ,r') } 


where   At   =    t-t1    and    Ar   =   r-r ' .       If   Eq.     (2.2-60)     is    substituted 
into  Eq.     (2.2-59),    then 


00  00 


// 

—  oo      —  x> 


R    (t,t,(r,r,l    =     I       I     Sx(f,v)RH(f,f,v>,v    ;    t,  t*,r ,  r')exp  (  +  j  2uf  At) 


exp(-j2TTvAr)dfdv/  (2.2-62) 


and  if  it  is  further  assumed  that  the  transfer  function  is  wide- 
sense  stationary  in  both  time  and  space,  then  Eq.  (2.2-62)  re- 
duces to 


oo  oo 


R    (At,Ar)    =      /       I     Sx(f ,v)RR(f ,f ,v,v    ;     At,Ar)« 

—  CO       —  oo 

exp(  +  j2Trf  At)exp(-j2irv'  Ar )  df dv_  (2.  2-63) 

which  implies  that  the  output  y(t,r)  is  also  wide-sense  sta- 
tionary in  both  time  and  space.   If  we  define  the  output  power 
spectral  density  S  (n,8_)  as 


S  (n,S)  =  F..F.  {R  (At,Ar) },       (2.2-64; 

y    —      At  Ar_  y      — 
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then   substituting   Eq.     C2.2-63)     into    the   R.H..S.    of   Eq.     (2.2-64) 
yields 


CD  CO 


// 

—  CO       -*oo 


S    (n,6.)    =     I      I     Sx(f ^vjRgtf ff ,v,v    ;    n-f/    8_-v)dfdv 

(2.2-65) 


where, from  Eq.     (2.2-41), 


00  oo 


// 

—  OO  —  OO 


RgCf^,^,^   ;    n-f/8_-v_)    =     I       J     RH(f/f/^v_   >    At,  Ar )  exp  [-j  2tt  ( n-f )  At] 


exp  [+j2tt  (S_-v_)  •  ArjdAtdAr. 

(2.2-66) 

Note  that  Eq.  (2.2-65)  is  in  the  form  of  a  multidimensional 
convolution  integral  analogous  to  Eq.  (2.1-46)  for  the  deter- 
ministic case.   The  convolution  process  accounts  for  the  fre- 
quency spreading  (frequency  in  HZ.  and  spatial  frequencies  in 
cycles/m. )  of  the  input  power  spectral  density. 

Example  2 . 2-1 

Using  Eqs.  (2.2-65)  and  (2.2-66),  let  us  compute  the 
output  power  spectral  density  for  the  case  when  the  transfer 
function  is  time- invariant,  space- invariant,  and  deterministic. 

First  note  that 

R„(f,f,v,v  ;  At,Ar)  =  E{H(f,v  ;  t,r)H*(f,v  ;  t',r')}. 
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If  H(_f,_v  ;  t,r)  is  hoth  time-invariant  and  space-invariant, 
then  [see  Eq.  (2.1-16)] 

H(f ,v  ;  t,r)  =  H(f,v) , 

and  if  the  transfer  function  is  also  deterministic,  then 


RR(f,f,v,v  ;  At,Ar)  =  E{ |H(f,v) |2} 

=  |H(f,v) |2. 


Therefore,    Eq.     (2.2-66)    becomes 


RgCf  ,f  ,v,v  ;  n-f  ,£-v)    =    |H(f,v)|2     /    exp[-j2ir  (n-f )  At]dAt 

—  CO 
00 

j     exp[+j2TT  (8_-v_)  •  ArjdAr 


=    |H(f,v) |25(n-f)6(S-v)/ 


and  substituting  this  expression  into  Eq.  (2.2-65)  yields  the 
desired  result 


S  (n,3)  =  |H(n,B) |2S  (n,3)  . 
y    —  —     x    — 

By  replacing  q  with  f  and  &  with  v_,    we  finally  obtain  [compare 
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with   Eq.     (2.1-48) J 


S     (f,v)    =    |H(f,v) |Sx(f,v)     . 


***** 
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III.    COUPLING  EQUATIONS 

We  will  now  present  those  equations  which  couple  the 
transmitted  and  received  electrical  signals  to  the  transfer 
function  of  the  ocean  medium  via  the  transmit  and  receive  far- 
field  directivity  functions.   Referring  to  Fig.  2,  the  fre- 
quency and  angular  spectra  of  the  input  acoustic  signal  to 
the  medium  and  the  output  electrical  signal  from  the  receive 
aperture  (array)  are  given,  respectively,  by  the  following 
equations  [2] : 


00 


xM(f'^)  ■   /  X(f,a)DT(f ,v-a)da  (3-1) 


and 


I   V"'i)D, 


Y(n,Y_)  =  |  Y„(n,3)DR(n,Y-e.)d6.  (3-2) 


where  the  frequency  and  angular  spectrum  of  the  output  acoustic 
signal  from  the  medium  is  given  by  [see  Eq .  (2.1-46)] 


CO      OO 


j  f  Vf^> 


YM("<A>  *  I   I  xM(f 'ji'BM(f'iL  ;  n-f,£-v)dfdv      (3-3) 


_  oo   —  oo 
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where    [see   Eq.     (2.1-32)] 


ao  oo 


BM 


(f,v_   ;    <}>,<_)    -     /       /     H    (f,v_   ;    t,r)  exp  C  — j  2ir«j>  t)  • 


•oo       —  oo 


exp(  +  j2TT<_«r)dtdr  (3-4) 

is  the  bi-frequency  function,  and  H  (f,v_  ;  t,r)  is  the  time- 
variant,  space- variant,  transfer  function  of  the  ocean  medium, 
Equations  (3-1)  thru  (3-4)  are  the  basic  coupling  equations. 
However,  w.r.t.  Eq.  (3-1),  we  have  the  following  additional 
relationships : 


oo 

fh 


XM(f,v)  =  FtFr{xM(t,r)  }  =   /  /  xM(t,r)exp(  +  j2Trv.r)dV. 


-oo  V 

exp(-j2Trft)dt  (3-5) 


where  [2] 


/ 


xM(t,r)  =   /  X(f  ,r)AT(f  ,r )  exp  (  + j  2irf  t)  df  ,         (3-6) 


t  r 

-oo"V 


X(f,a)  =  F+.Fy.{x(t,r)  }  =  /  /   x  ( t ,  r )  exp  (  +  j  2™  -r)  dV- 

exp(-j2irft)dt  (3-7) 
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is  the  frequency  and  angular  spectrum  of  the  transmitted  elec- 
trical signal,  and 

DT(f,a)  =  Pr{AT(f/r)}  =  /  AT(f,r)exp(  +  j2TTa-r)dV 

V 

(3-8) 

is  the  far-field  directivity  function  (beam  pattern)  of  the 
complex  transmit  aperture  A  (f,r).   Also,  w.r.t.  Eq.  (3-6), 


oo 

■/ 


X(f,r)    =  Ft(x(t,r)}   =     /    x(t,r)exp(-j2irft)dt 


3-9) 


and 


AT(f,r)    =    F^X{DT(f,a)}    =      /     DT  (f  ,  a)  exp  (-  j  2tto_  -r )  da    . 

—  oo 

(3-10) 


Similarly,  w.r.t.  Eq.  (3-2),  we  have  the  following 
additional  relationships: 


Y(n,Y)  =  FfcF 


CO 

r(y(t,r)}  =  /  /  y(t,r)exp(  +  j2iry-r)dV- 


exp(-j2irnt)dt  (3-11) 
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where    [2] 


f    V 


y(t,r)    =     /    Y    Cn,r)A    (n,r)expC+j27rnt)dn/  (3-12) 


the    frequency   and   angular    spectrum   of    the   output   acoustic    sig- 
nal   from   the   medium   is   given  by 


00 

YM(n,S)    =   FtFr{yM(t,r)  }    =     f  f  yM(tfr)  exp  (  +  j  2ir£-  r)  dV- 

—  oo      V 

exp(-j2Trnt)dt  (3-13) 


where    [see   Eq. (2.1-25)] 


00  00 


yM(t,r)    =      I       I    XM(f,v)HM(ffv    ;    t,r)  exp  (  +  j2irf  t) 


—  oo       — oo 


exp(-j2TTvr)dfdv    ,  (3-14) 


and 


■/ 


DR(i'i)    =   Fr{AR(n/£)}    =    /   AR(n,r)exp(  +  j27T8_-r)dV 


(3-15) 


is  the  far-field  directivity  function  (beam  pattern)  of  the 
complex  receive  aperture  A _(n#r).   Finally,  w.r.t.  Eq.  (3-12), 


58 


CO 

/ 


YM(n,r)    =  Ft{yM(t,r)}    =    |     yM(t,r)  exp  (-j27rn  t)  dt 

(3-16) 


and 


CO 

/v- 


AR(n  »r)    =  F~  {DR(n  ,g_)}  = 


(3-17) 
Substituting  Eq.  (3-14)  into  Eq.  (3-16)  yields 


CO  CO 


YM(n'£)    =     I       I     xM(f '^)BM(f '-  ;    n"f >£)exP(-J2™-r)dfdv_ 


—  CCj  —CO 


(3-18) 

where 

oo 

BM(f,\i   ;    <j>,r)    =     /    HM(f,v    ;    t,r)  exp  (-j  2tt<J)  t)  dt 

—  00 

(3-19) 

which  is  not  the  same  as  the  bi-frequency  function  [see  Eq. 
(3-4)] . 

The  components  of  the  vectors  <* ,  v_,x  6_,  and  y_   are  spatial 
frequencies  with  units  of  cycles/m.  and  r  =  (x,y,z). 

A  moments  reflection  leads  one  to  the  conclusion  that 
Eq.  (3-12)  is  well  suited  for  space-time  signal  processing 
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applications.   By  inspecting  Eqs.  C3-18)  ,  C3-19)  ,  and  (.3-1)  , 
it  can  be  seen  that  the  output  signal  given  by  Eq.  (.3-12)  can 
be  expressed  as  a  function  of  time  and  space  in  terms  of  the 
transmitted  electrical  signal,  the  transmit  aperture,  the 
transfer  function  of  the  ocean  medium,  and  the  receive  aperture 
Equation  (3-2),  however,  represents  the  theoretical  output 
"spectrum",  albeit  dependent  upon  the  same  system  functions. 

If  the  ocean  medium  is  modelled  as  a  random  filter, 
then  y(t,r)  given  by  Eq.  (3-12)  is  a  random  process  with  auto- 
correlation function 


Ry(t,f,r,r')  =  E{y(t,r)y*(t,,r')  }.        (3-20) 


It  is  probably  obvious  by  now  that  the  weak  link  in  the 
coupling  equations  is  the  ocean  medium  transfer  function 
H  (f,v_  ;  t,r)  and  its  corresponding  autocorrelation  function 
[see  Eq.  (2. 2-2) ] 


Rjj  (f,f',v,v'  ;  t,t\r,r')  =E{HM<f,v  ;  t,r)H*(f',v'  ;  f,^)} 
M 

(3-21) 


If  one  is  not  able  to  specify  a  realistic  functional  form  for 
H  (f,v_  ;  t,r) ,  then  the  equations  presented  in  this  section 
represent  an  exercise  in  the  methods  of  linear  system  theory. 
This  leads  us  to  the  main  purpose  of  this  paper  which  is  to 
demonstrate  the  derivation  of  an  ocean  transfer  function  based 
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upon  the  W.K.B.  approximation.   The  derivation  is  discussed  in 
Section  4.1. 

Finally,  note  that  the  autocorrelation  function  given 
by  Eq.  (3-21)  can  be  referred  to  as  a  generalized  coherence 
function  since  it  is  a  generalization  of  the  two-frequency  cor- 
relation function  or  two-frequency  mutual  coherence  function 
based  upon  linear,  time-variant  filter  theory  as  discussed  by 
Ishimaru  [51],  for  example.   Equation  (3-21)  provides  informa- 
tion concerning  the  amount  of  spreading  in  time-delay,  space, 
frequency,  and  spatial  frequencies  that  a  transmitted  signal 
will  be  expected  to  experience  as  it  propagates  in  a  random 
medium.   The  result  of  all  this  spreading  is,  of  course,  dis- 
tortion in  pulse  shape. 
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IV.     ANALYSIS 

4 . 1   Transfer  Function 

We  will  now  proceed  to  derive  a  transfer  function 
which  models  the  bistatic  communication  channel  geometry  shown 
in  Fig.  1.   The  communication  channel  is  regarded  to  be  the 
ocean  volume  between  the  apertures  so  that  surface  and  bottom 
scattering  effects  are  not  included.   Both  apertures  are  sta- 
tionary (not  in  motion) ,  and  it  is  assumed  that  no  discrete 
point  scatterers  (such  as  bubbles,  fish,  etc.)  are  in  the 
volume  between  the  apertures.   No  motion  implies  that  the  re- 
sulting transfer  function  will  be  time- invariant. 

The  propagation  of  small  amplitude  acoustic  signals 
in  the  ocean  from  the  transmit  aperture  to  the  receive  aperture 
can  be  described  by  the  following  linear,  inhomogeneous ,  scalar 
wave  equation: 


V2*(t,r) 2^ —  "^2  *(t'£>  =  XM(t'-} 

c  (r)   3t 

(4.1-1) 


where  <P(t,r)  is  the  velocity  potential  at  time  t  and  position 
r  =  (x,y,z),  x„(t,r)  is  the  source  distribution  [see  Fig.  2 
and  Eq.  (3-6)],  and  c (r )  is  the  speed  of  sound  in  the  ocean. 
Since  the  coupling  equations  discussed  in  Section  III  already 
allow  for  an  arbitrary  x  (t,r)  with  corresponding  frequency 
and  angular  spectrum  SL-Cf , v) ,  we  need  only  find  the  solution 
to  the  following  Helmholtz  wave  equation: 
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V2*tr)    +   k2   n2(r)*(.r)    =    0 


where 


C4.1-2) 


k      =    27Tf/c      =    2tt/a 
o  o  /    o 


(4.1-3) 


is  the  constant,  reference  wavenumber, 


n(r)  =  cQ/c(r) 


(4.1-4) 


is  the  random  index  of  refraction, 


c   =  c(r  )  =  fX 
o     — o      o 


(4.1-5) 


is  the  constant,  reference  speed  of  sound  at  the  source  posi- 
tion r   =  (x  ,y  ,  z  ) ,  and 
— o     o  J  o     o 


*(t,r)  =  *(r)  exp(  +  j2TTft) 


(4.1-6) 


is  the  time-harmonic  solution  of  Eq.  (4.1-1)  when  x  (t,r)  is 
set  equal  to  zero, and  where  ^(r)  is  the  solution  of  Eq.  (4.1-2). 
Note  that  the  wavenumber 


k(r)  =  2irf/c(r) 


(4.1-7) 


can  be  expressed  as 


k(r)  =  kQn(r) 


(4.1-3) 
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Therefore.  k(r  )  =  k   since  n(r  )  =  1.   The  index  of  refraction 
— o     o         — o 


is  commonly  written  as  [40,  42-45] 


n(r)  -  nD(r)  +  nR(r)  (4.1-9) 

or 

n(r)  =  nD(r)  +  a(r)nNR(r)  (4.1-10) 


where  n_ (r )  is  the  deterministic  component  and  is  usually 
close  to  unity  in  value,  n_(r)  is  the  random,  zero  mean  compo- 
nent,  a (r )  is  the  standard  deviation  of  nR(r) ,  and 


nNR(r)  =  nR(r)/a(r)  (4.1-11) 


is  the  normalized  random  component  with  zero  mean  and  variance 
equal  to  unity.  We  shall  work  with  Eq.  (4.1-10)  in  this  paper 
Note  that  the  average  value  of  n(r_)  is  equal  to  n_(r). 

Let  us  assume  that  the  speed  of  sound  is  only  a 
function  of  the  depth  y,  i.e.,  c (r )  =  c(y),  so  that  Eq.  (4.1-21 
reduces  to 


72*(x,y,z)  +  k^  n2(y)*>(x,y,z)  =  0      (4.1-12) 


where,  from  Eq.  (4.1-10), 


n(y)  =  nD(y)  +  a(y)nNR(y).  (4.1-13) 
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By  using  the  method  of  separation  of  variables  and  the  W.K.B. 
approximation  [46,47],  an  approximate  solution  of  Eq.  (4.1-12 
is  given  by 

*>(x,y,z)  »  exp[-jkx(x-xQ)]  a^.  (y) exp[  +  je y(y) ] exp[-jkz (z-zq) J 

(4.1-14) 


where 


aY(y)   =  [kY(y)]"1/2         (4.1-15) 


and 


8Y(y)  =  expj-j  /   kY(?)d?[  (4.1-16) 

yo 


1/2 
ky(y)  =  [k2n2(y)  -  k^-k2]  (4.1-17) 


where  k  ,  k  (y)  ,  and  k   are  the  components  of  the  propagation 
X    i  ^ 

vector 


k(y)  =  kxx  +  ky(y)y  +  kzz.        (4.1-18) 

Note  that  kx  and  k   are  constants  while  k„ (y)  is  a  function  of 

the  depth  y.   In  addition,  Eqs.  (4.1-14)  thru  (4.1-16)  allow 

for  a  general  source  location  with  position  vector  r   =  (x  ,y  ,  z  I 
J  c  — o     o   o  o 

The  W.K.B.  approximation  given  by  Eqs.  (4.1-15)  and 
(4.1-16)  is  a  valid  solution  for  the  depth  dependence  provided 


65 


that  [46,47,50]  1)  the  transmitted  frequencies  f  are  high,  2) 
the  sound  speed  profile  c(y)  is  slowly  changing,  and  3)  the 
depth  interval  from  y   to  y  does  not  include  any  turning 
points.   A  turning  point  exists  at  y  =  yT  if  k  (y_)  =  0. 


Since  the  reference  propagation  vector  k   can  be 


expressed  as 


then 


where 


and 


where 


and 


k^  =  kxx  +  kyy  +  kzz  ,  (4.1-19) 


ko  =  1 2So  I  2  =  kX  +  kY  +  kZ        (4.1-20) 


kx  =  koUQ  (4.1-21) 


kv  =  k  v  (4.1-22) 

Y     O  O 


k7  =  k  w  (4.1-23) 

z    o  o 


u   =  sine  ccsiij  (4.1-24) 

o       o    ro 


v  =  sine  sinib   =  coss  (4.1-25) 

o       o    o       o 


w   =  cose  (4.1-26) 

o       o 


are  the  direction  cosines  w.r.t.  the  positive  X,Y,  and  Z  axes, 

respectively,  and  (e  ,ti>    )  are  the  vertical  and  azimuthal 

o   o 

spherical  angles  measured  w.r.t.  the  positive  Z  and  X  axes, 
respectively,  representing  the  initial  directions  of  wave  pro- 
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pagation  (see  Fig.  7).   From  Eq.  (4.1-20)  we  can  write  that 


kX  +  kZ  =  ko  "  kY  (4.1-27) 


and  substituting  Eqs.  (4.1-27)  and  (4.1-22)  into  Eq.  (4.1-17) 
yields 

1/2 


kY(y)  =  K,    jl  +  fn2(Y^n 


v 
o 


(4.1-28) 


Note,  that  k„ (y  )  =  k   since  n(y  )  =  1.   Therefore,  the  square 
root  expression  in  Eq.  (4.1-28)  is  responsible  for  changing 
the  initial  direction  of  propagation  ky.   If 


[n2(y)  -  l]/v^|  <  1  ,  (4.1-29) 


then  the  square  root  expression  in  Eq.  (4.1-28)  can  be  approx- 
imated by  the  first  two  terms  in  a  binomial  expansion  yielding 


ky(y)  «  kY  +  k2[n2(y)  -  l]/(2ky)  .        (4.1-30) 


It  will  be  shown  later  that  the  binomial  expansion  criterion 
given  by  Eq.  (4.1-29)  can  be  related  to  the  critical  angle  of 
incidence,  and  hence,  total  reflection.   The  assumption  con- 
cerning the  absence  of  turning  points  will  also  be  discussed 
later,  but  first  let  us  return  to  the  derivation  of  the  trans- 
fer function. 
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Recall  that  we  designated  the  vector  v_  to  represent 
the  transmitted  spatial  frequencies  as  in  X  (f,v_)  (see  Fig.  2). 
If  we  let 


then 


and 


v  =  (fx,fy,fz) ,  (4.1-31) 


fv  =  kv/(27r)  =  u/Xrt  (4.1-32) 

A       A  O    O 


fy  =  ky/(2TT)  =  vQ/Xo  (4.1-33) 


fz  =  kz/(2Tr)  =  wQ/Xo  (4.1-34) 


are  the  spatial  frequencies  in  the  X,Y,  and  Z  directions,  re- 
spectively, where  c   =  c(y  )  =  fX  .   Therefore,  if  Eqs.  (4.1-30) 
thru  (4.1-34)  are  substituted  into  Eqs.  (4.1-14)  thru  (4.1-16), 
then  Eq.  (4.1-6)  becomes 


-j[k£/(4»£Y)]  /   [n2 


I* 


*(t,x,y,z)  « 


yo 


2    2  | 1/2 

2irfy  +  k^  [nZ(y)  -  1]  /(4^fy) 


exp(+j27rft)exp  [-j2ttv_-  (r-r  )].      (4.1-35) 


Since  the  input  to  the  communication  channel  is  the  time-har- 
monic plane  wave 


exp  (  + j2iTf  t)  exp  [- j2t:v_'  (r-r  )], 
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and  since  Eq.  (4.1-35)  represents  the  output  from  the  channel 
at  any  time  t  and  position  r,  then  [see  Eq.  (2.1-21)] 

HM(f,v  ;  t,r)  =  HM(f,v  ;  r)  =  HM(f,fY  ;  y) 

(4.1-36) 

where  the  random,  time- invariant,  space-variant  transfer  func- 
tion of  the  ocean  medium  is  given  by 


HM(f,fy  ;  y)  =  -Vf,fY  ;  y)exp[+jOM(f,fY  ;  y)  ] 


where 


(4.1-37) 


Vf'fY  ;  y)  =  {27TfY  +  ^[n2(y)-l]/(41rfY)}"1/2 


(4.1-38) 


6M(f,fY  ;  y)  =  -[k^/(4TrfY)]  /   [n*(c)-l]dC 


(4.1-39) 


k   =  2-rrf/c 
o        o 


(4.1-3) 


and  from  Eq.  (4.1-13), 


2        2  2     2 

n  (y)  =  nDCy)  +  2nQ (y) a (y) nNR(y)  +  a  (y)nNR(y) 


(4.1-40) 


If  the  medium  is  "weakly  irregular"  or  "weakly  scattering" , 
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i.e.,  if  a  (y)  is  very  small  compared  to  unity  so  that 


ko°  Cy)sY  <<  1  (4.1-41) 


where  s   is  the  scale  size  of  the  irregularities  (i.e.,  the 

average  distance  over  which  the  refractive  index  fluctuations 

2 

remain  correlated  [48]),  then  terms  involving  a    (y)  can  be  ne- 
glected [42]  and  Eq.  (4.1-40)  reduces  to 


n2(y)  M  n2(y)  +  2nD (y) a (y) nNR(y) .   (4.1-42) 


Therefore,  by  using  Eq.  (4.1-29)  to  simplify  Eq.  (4.1-38)  and 
upon  substituting  Eq.  (4.1-42)  into  Eq .  (4.1-39),  one  obtains 
the  following: 


and 


Vf'fY  ;  y)  =  (27rfY)'1/2  (4.1-43) 


eM(f,fY  ;  y)  =  eMD(f,fY  ;  y)  +  °MR(f'fY  ;  y) 


(4.1-44) 
where  y 

0MD(f,fY  ;  y)  =  -[^/(4TTfY)]  J       [n2(c)-l]d?     (4.1-45) 

Yo 

is  the  deterministic  or  average  component  of  the  phase  function, 
and 
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0MRCf,fY  ;  y)  =  -[k^/(2TTfy)]  J      nDU)aU)n^RU)dz 


yo  (4.1-46) 

is  the  random  component. 

The  medium  transfer  function  H   given  by  Eq.  (4.1-37) 
with  amplitude  and  phase  functions  A^  and  0   specified  by  Eq. 
(4.1-43)  and  Eqs.  (4.1-44)  thru  (4.1-46),  respectively,  in- 
dicate that  for  a  weakly  irregular  medium,  the  major  effect  of 
the  medium  is  to  angle  modulate  the  transmitted  field.   Further- 
more, if  nD(y)  a  1,  then  0   (f,f   ;  y)  «  0   [see  Eq.  (4.1-45)] 
and,  as  a  result,  the  angle  modulation  is  due  strictly  to  the 
random  fluctuations  <j(y)nN_(y)  of  the  index  of  refraction  [see 
Eq.  (4.1-46)].   Note  that  the  angle  modulation  process  is  often 
referred  to  as  "scattering"  [48].   Also  note  that  the  transfer 
function  derived  in  this  section  can  be  written  as  the  product 
of  two  functions,  one  deterministic  and  the  other  random,  i.e., 

HM   =     [VXp(  +  jW][eXp(  +  j°MR)] 

which  agrees  v/ith  the  assumed  transfer  function  expressions  of 
Laval  and  Labasque  [36]  and  with  the  general  practice  of  re- 
presenting a  field  in  a  random  medium  as  the  product  of  a  de- 
terministic and  a  random  function  [40] . 

Let  us  now  discuss  the  two  important  assumptions 
responsible  for  the  derivation  of  the  transfer  function,  namely, 
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the  aforementioned  binomial  expansion  criterion  given  by  Eq. 
(4..  1-29)  and  the  assumption  regarding  the  absence  of  turning 
points. 

Since 

n(y)  =  cQ/c(y) ,  (4.1-47) 

substituting  Eqs.  (4.1-25)  and  (4.1-47)  into  Eq.  (4.1-29) 
yields 


(|c2(y)  -  c2l)1/2 

DS    J    '     J    O  | 

c(y) 


(4.1-48) 


where  $   is  the  angle  of  incidence  of  the  reference  propagation 

vector  k   (see  Fig.  7) .   Now  recall  that  the  critical  angle  of 

incidence  s   associated  with  a  time -harmonic  plane  wave  inci- 
c 

dent  upon  a  plane  boundary  between  two  fluid  media  is  given  by 
[49]  (see  Fig.  8) 


sin6c  =  Cl/c2 


c2  >  c 


(4.1-49) 


or,  equivalently, 


=  cos 


-1 


i    2  2\ 

(c2  -  C]_) 


1/2 


c2  >  c± 


(4.1-50) 


where  c~  must  be  greater  than  c,  for  3   to  exist.   When  the 

angle  of  incidence  B  >  3   /  there  is  total  reflection  and  hence, 

c  
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Fig.  7.   The  reference  propagation  vector 
k   and  associated  angles. 
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Fig.  8.   Illustration  of  the  angle  of  incidence 

B  and  the  critical  angle  of  incidence  Sc 
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no  transmission  of  energy  into  medium  II.   The  speed  of  sound 
c (y)  generally  increases  with  depth,  for  example,  in  the  deep 
ocean.   Also,  c(y)  increases  above  and  below  the  SOFAR  channel 
axis.   Therefore,  when  c(y)  >  c  ,  the  absolute  value  sign  in 
Eq.  (4.1-48)  can  be  removed,  and  by  comparing  Eq.  (4.1-50) 
with  the  R.H.S.  of  Eq.  (4.1-48),  the  binomial  expansion  crite- 
rion indicates  that  the  initial  angles  of  transmission  3   must 

o 

be  less  than  the  "critical  angle"  in  order  to  avoid  total  re- 
flection, and  thus,  passing  thru  a  turning  point  [47,50]. 
This  is  consistent  with  the  fact  that  the  W.K.B.  approximation 
is  invalid  at  a  turning  point  [47,50] . 
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4 . 2   Output  Electrical  Signal 

Now  that  we  have  derived  an  ocean  medium  transfer 
function,  let  us  demonstrate  the  use  of  the  coupling  equations 
presented  in  Section  III  by  calculating  the  output  electrical 
signal  y(t,r)  from  the  receive  aperture  (array)  as  given  by 
Eq.  (3-12) .   As  can  be  seen  by  inspecting  the  integrand  of 
Eq.  (3-12),  we  need  expressions  for  both  of  the  kernels  Y„(n,r) 
and  A_(n,r) .   Let  us  compute  YM(n/r)  first  [see  Eq.  (3-18)]. 

Assume  that  the  transmit  aperture  depicted  in  Fig.l 
is  a  planar  array  of  MxN  (odd)  complex  weighted  point  sources, 
centered  at  (x  ,y  ,  z  )  and  parallel  to  the  XY  plane.   In  addi- 
tion,  assume  that  the  complex  weights  are  separable.   Since  in 
most  practical  situations  an  identical  input  electrical  signal 
is  applied  to  all  elements  in  the  transmit  array  before  the 
complex  weights,  i.e.,  since 

x(t,r)  e  x(t) , 

then  Eq.  (3-7)  reduces  to 

X(f,a)  =  X(f) 5(a)  (4.2-1) 

since  F  {1}  =  5(a).   Substituting  Eq.  (4.2-1)  into  Eq.  (3-1) 
yields 

XM(f,v)  =  X(f)DT(f,v)  (4.2-2) 
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where 


(M-D/2  (N-D/2 

DT(f,v)    =  2md  2^1  cmdnexp(+j27rfxmdx)  • 

m=- (M-D/2   n=- (N-D/2 

exp  (+ j  2:rf Yndy)  exp  (  + j  27rf xx   ) 
exp  (  +  j2ufYyo)  exp  (+j  2irf  zzq) 

(4.2-3) 


is  the  far-field  beam  pattern  of  the  transmit  array,  c   and  d 

r  J        m      n 

are  complex  weights,  d   and  dY  are  the  interelement  spacings  in 
the  X  and  Y  directions,  respectively,  and  the  last  three  expo- 
nentials are  phase  factors  which  account  for  the  array  being 
centered  at  (x  , y  ,z  )  instead  of  at  the  origin  (0,0,0). 

Since  the  transfer  function  derived  in  Section  4.1 
is  time-invariant,  Eq.  (3-19)  becomes 


BM(f,v  ;  <J),r)  =  HM(f,v  ;  r)5(cj>)        (4.2-4) 


and,  as  a  result, 


BM(f'^  J  n-ff.r)  =  HM(f,v  ;  r)6(n-f).       (4.2-5) 


Therefore,  substituting  Eqs.  (4.2-2)  and  (4.2-5) 
into  Eq.  (3-18)  yields 
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YM(n'E>  =  X(n)  /  DT(n  ,v_)HM(n,v_  ;  r)exp(-j2*v/ r)dv_, 

i-  go 

(4.2-6) 
and  upon  substituting  Eq.  (4.2-6)  into  Eq.  (3-12),  we  obtain 


/ 


y(t,r)  =  /  X(f)   /  D_(f,v)H  (f,v_  ;  r )  exp  (-  j2^v_-  r )  dv/ 


M 


A_(f  ,r)exp(  +  j27Tft)df  (4.2-7) 

is  — 


where  n  was  replaced  by  f  since  there  is  no  frequency  (HZ) 
spreading  in  this  example  due  to  the  time- invariant  property 
of  the  transfer  function.   We  need  to  specify  the  complex 
aperture  function  A_(f,r)  next. 

Assume  that  the  receive  aperture  depicted  in  Fig.  1 
is  a  planar  array  of  M'xN1   (odd)  complex  weighted  point 
sources,  centered  at  (x^y-,,  z_)  and  parallel  to  the  XY  plane. 

i\   i\   K. 

In  addition,  assume  that  the  complex  weights  are  separable. 
Therefore,  the  receive  aperture  function  is  given  by 

(M'-l)/2     (N'-l)/2 
AR(f,r)  =      X)  Z     c!d'5(x-[xR  +  id']>- 

i=-(M'-l)/2  q=-(N'-l)/2 

My-[yR  +  qd£])  6(z-zR) 

(4.2-8) 
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where  c!  and  d1  are  complex  weights  and  d'  and  d'  are  the  inter- 
1      q        c  X      Y 

element  spacings  in  the  X  and  Y  directions,  respectively. 

Upon  substituting  Eqs.  (4.1-31),  (4.1-36),  (4.2-3), 
and  (4.2-8)  into  Eq,  (4.2-7)  and  recalling  that  r  =  (x,y,z), 
one  obtains 


(M'-l)/2     (N'-l)/2 

00 

y(t,x,y,z)  =      ^  ^2  cidq  /  X(f) 

i=-(M'-l)/2  q=-(N'-l)/2 


(M-D/2 


Z    Cm  (    exP(-^fXAXim)dfx' 


m=- (M-D/2 


(N-D/2 


^2 


where 


and 


n=- (N-D/2     -co 

CO 

/  exp(-j2iTf  AZ)df  exp(  +  j2Trft)df • 

—  CO 

5(x-[xR+  id^])5(y-[yR+  qd^])5(z-zR)     (4.2-9) 


AX.   =  (x  -x  )  +  (id'  -  mdv)         (4.2-10) 


AYqn  =  (yR'Yo)  +  (qdY  "  ndY}         (4.2-11) 


AZ  =  z„  -  z  .  (4.2-12) 

R     o 
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If  we  now  change  variables  from  spatial  frequencies  to  direc- 
tion cosines  by  substituting  Eqs.  (4.1-32)  thru  (4.1-34)  into 
Eq.  (4.2-9)/  and  treat  the  frequency  variable  f  as  a  constant 
with  respect  to  the  spatial  frequency  integrations,  then  Eq. 
(4.2-9)  becomes 

(M'-l)/2     (N'-l)/2 


D  E      Cidq  / 


y(t,x,y,z)  -  (2/c^)      zl  5"!    c^dA  /  f3*< 


i=-(M'-l)/2  q=-(N'-l)/2 


(M-D/2 


Zc  si 
m 

m=- (M-D/2 


nc(2fAX.  /c  )i 

m         lm  o  ' 


(N-D/2      +1 


E    dn  /  HM(f'  tfvo/cQ];  yR+  qdj) 


n=-(N-l)/2     a 

q 


exp(-j2Tr[fv   /c    ]  AY      )  dv 
r      J       L      o      oJ      qn        c 


(b   exp  (-jirb    f  AZ/c    )sinc(b   fAZ/c    )} 
q  q  o  q  o 


exp(  +  j2Tift)df  5(x-     [xR+    id£])5(y-     [yR+   qd^])5(z-zR) 


(4.2-13) 


since  in  our  example  problem 


and 


where 
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-1  <    u     <   +1  (4.2-14) 

o 


a    <    v     <  +1  C4.2-15) 

q  o 


0  <  wQ  <    b  (4.2-16) 


aq    =     (lnD(yR   +    qdY}     "    1I)V2  (4.2-17) 


bq   =    (1    -    |n^(yR  +   qd^)    -    l|)1/2       (4.2-18) 


smc(x)    = (4.2-19) 

ux 


and  n  (•)  was  replaced  by  n2 ( • )  in  Eqs.  (4.2-17)  and  (4.2-18). 
The  transfer  function  H   is  given  by  Eqs.  (4.1-37)  and  (4.1-43) 
thru  (4.1-4  6)  with  the  exception  that  the  lower  limit  of  in- 
tegration y   in  Eqs.  (4.1-45)  and  (4.1-46)  must  be  replaced  by 
y  +  nd  which  is  the  Y  coordinate  of  a  point  source  in  the 
transmit  planar  array.   Note  that  if  nn ( • )  <  1,  then 


bq  =  nD(yR+  qd£).  (4.2-20) 


Equation  (4.2-13)  is  the  desired  result,  i.e.,  it  represents 
the  random  output  electrical  signal  from  each  element  in  the 
receive  array  in  terms  of  the  transmitted  electrical  signal, 
the  transmit  and  receive  arrays,  and  the  random  transfer  func- 
tion of  the  ocean  medium.   If  y(t,x,y,z)  given  by  Eq.  (4.2-13) 
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is  complex,  then  simply  take  the  real  part  to  obtain  the  real 
output  electrical  signal. 
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4 . 3   Coherence  Function 

Upon  inspecting  Eq.  (4,2-13),  it  can  be  seen  that 
the  autocorrelation  function  of  the  output  electrical  signal 
will  depend  upon  the  autocorrelation  function  of  the  transfer 
function 


Rjj  (f,f»,fY,f£  ;  y,y')  =  E{HM(f,fy  ;  y)H*(f\fy  ;  y»)) 

M 

(4.3-1) 


which  is  also  known  as  the  coherence  function.   Substituting 
Eqs.  (4.1-37),  (4.1-39),  and  (4.1-43)  into  Eq.  (4.3-1)  yields 

RR    (f,f',fy,fy    ;    y,y')    =    [271  (fyfy)  1/2J         E   exp  (  +  j  [KI  (y)  +K  •  I  (y  ' )  ]  ) 


•         •     *  -     X  -  -  -     -  L  x    X  J  1 

M  I 


(4.3-2) 


where 


K  =   -   k^/(4Trfy)  (4.3-3) 


K'    =      +(k^)2/(47tfy)  (4.3-4) 


y 

[n2(         -    l]d?  4.3-5 


.    * 


y, 


y' 

I(y')    =  [n2(r)    -    i]dc  (4.3-6) 


')    -     f        [n: 


yo 
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k   =  2Trf/c  (4.3-7) 

o        o 


k'  =  2Trf/c  (4.3-8) 

o         o 


2 

and  n  (•)  is  given  by  Eq.  (4.1-4  2)  where  the  expectation  appear- 
ing in  Eq.  (4.3-2)  is  the  characteristic  function  of  the  random 
quantity  [KI (y)  +  K'I(y')].   If  it  is  assumed  that  I (y)  is  a 
real  Gaussian  random  process,  which  implies  that  the  index  of 
refraction  is  a  real  Gaussian  random  process,  then  Eq.  (4.3-2) 
can  be  written  as 


Kg  (f,f',f  ,f£  ;  y,y')  =  [2*  |f  f«]  1/2J  exp(+j  [^(f  ,fy  ;  y)- 

M 

0MD(f,f'  ;  y')])- 

exp(-E{0^R(f,fY  ;  y)  }/2)  • 

exp(+E{9MR(f,fY;  y)3MR(f',f'  ;  y')» 

exp(-E{0^R(f,fY  ;  y')}/2)    (4.3-9) 


where 
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and 


(4.3-12) 


R         (y,y')    =   Edi^yJn^y')}.         (4.3-13) 


If  it  is  further  assumed  that  the  deterministic  component  n_(y) 
of  the  index  of  refraction  is  equal  to  unity  [40,  42-45],  and 
that  the  random  component  n_ (y) ,  and  hence,  the  normalized 
random  component  nN_(y)  is  wide-sense  stationary,  i.e., 


R   (y.y1)  =  R   (Ay) 

^TR  NR 


(4.3-14) 


where    Ay  =   y-y',    then    [see   Eq.  (4.1-45)] 


0MD(f'fY  '  y>  -  0 


(4.3-15) 
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and 


0MD(f''fY  '*  y'1  =  °' 


(4.3-16) 


and  Eqs.  (4.3-10)  thru  (4.3-12)  become,  respectively, 
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where  a  is  the  constant  standard  deviation  of  the  Gaussian, 
zero  mean,  wide-sense  stationary,  random  component  n_  (y)  of 
the  index  of  refraction. 
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V.      SUMMARY  and  CONCLUSIONS 

A  consistent  notation,  fundamental  input-output  relations, 
and  various  time-space  transformations  for  both  deterministic 
and  random  linear,  time-variant,  space-variant  filters  have 
been  established.   The  notation  is  consistent  in  the  sense  that 
all  of  the  various  input-output  relations  which  are  based  upon 
the  general  theory  will  reduce  to  the  classical  relations  of 
linear,  time- invariant  filter  theory.   These  results  should  be 
of  interest  to  persons  involved  in  the  general  area  of  linear 
systems  theory,  and  not  only  to  those  involved  in  underwater 
acoustics. 

With  the  use  of  the  method  of  separation  of  variables 
and  the  W.K.B.  approximation,  a  mathematical  expression  of  a 
time-invariant,  space-variant,  random  transfer  function  of  the 
ocean  medium  was  derived.   The  transfer  function  was  time-in- 
variant instead  of  time-variant  because  motion  was  not  consi- 
dered in  the  present  derivation.   The  transfer  function  modelled 
the  ocean  volume  between  transmit  and  receive  apertures  (arrays) . 
The  ocean  volume  was  characterized  by  a  random  index  of  refrac- 
tion (sound  speed  profile)  which  was  a  function  of  depth.   The 
index  of  refraction  was  decomposed  into  deterministic  and  random 
components. 

In  addition  to  the  transfer  function  derivation,  two 
example  calculations  were  made.   The  first  example  demonstrated 
the  use  of  the  coupling  equations  and  involved  the  derivation 
of  a  mathematical  expression  for  the  random  output  electrical 
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signal  at  each  element  in  a  receive  planar  array  of  complex 
weighted  point  sources.   The  output  signals  were  expressed 
in  terms  of  the  frequency  spectrum  of  the  transmitted  electrical 
signal,  the  transmit  and  receive  arrays,  and  the  previously 
derived  transfer  function  of  the  ocean  medium.   The  first  exam- 


ple demonstrated  that  an  output  electrical  signal  could  be 
derived  in  a  logical  and  straightforward  fashion.   The  second 
example  involved  the  derivation  of  the  coherence  function,  i.e., 
the  autocorrelation  function  of  the  transfer  function.   In 
order  to  obtain  somewhat  simplified  results,  it  was  necessary 
to  assume  that  the  random  component  of  the  index  of  refraction 
was  Gaussian  and  wide-sense  stationary. 
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